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Summary: 

This  final  technical  report  summarizes  the  progress  achieved  during  the 
past  year,  which  added  to  the  work  of  the  previous  four  years.  Progress 
achieved  during  the  first  four  years  of  the  project  were  summarized  in  the  Final 
Technical  Report  for  Contract  #  N00014-86-K-0791  (to  The  University  of 
Miami). 

The  primary  effort  during  the  past  year  focussed  on  analysis  of  data 
obtained  during  a  field  experiment  conducted  at  Cape  Canaveral  during 
Summer,  1987.  Analysis  was  pursued  using  a  new  approach  based  upon 
examining  the  spatial  and  temporal  characteristics  of  individual  sand  suspension 
events,  as  well  as  their  statistical  characteristics.  The  most  interesting 
conclusions  fell  along  two  lines  of  investigation:  1 )  the  presence  of  wave 
groups  enhances  the  suspension  of  sand  above  that  which  would  occur  with 
randomly  distributed  waves  with  the  same  height  distribution;  and  2)  the 
primary  difference  between  suspension  events  which  cause  erosion  and  those 
which  cause  accretion  relates  to  the  phase  of  the  wave  in  which  the  sediment 
is  suspended  off  the  seabed. 

Publications  describing  these  and  other  findings  are  listed  below,  and 
refereed  publications  are  included  in  the  appendix. 
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1.  Introduction 

The  movement  of  sediment  by  waves  in  coastal  regions  has  long  been  considered  to 
be  a  process  of  significance  to  short  term  bathymetric  changes  and  to  long  term 
coastal  erosion.  For  example,  there  is  a  prevalent  seasonal  trend  for  local  storm 
waves  to  erode  beaches  and  for  swell  waves  to  rebuild  beaches.  However,  the 
detailed  processes  by  which  waves  move  sand  are  not  well  known,  and  quantitative 
direct  measurements  of  moving  sediment  are  difficult  to  obtain.  In  general,  a  portion 
of  the  sediment  moves  in  suspension,  with  its  weight  supported  by  the  action  of 
fluid  turbulence.  Aspects  of  the  transport  of  suspended  sediment  are  treated 
elsewhere  in  this  volume  by  Cacchione  and  Drake  (Chapter  -21),  Kraus  and 
Horikawa  (Chapter  22).  and  Mehta  and  Dyer  (Chapter  23). 

Recent  observ’ations  indicate  that  near-bed  intermittent  suspension  is  an  identi¬ 
fiable  and  significant  process.  This  chapter  will  describe  the  structure  of  such  events 
observed  in  the  field  through  the  use  of  an  acoustic  profiling  sensor.  The  phe¬ 
nomenon  of  intermittent  suspension  has  been  observed  by  several  investigators 
including  Brenninkmeyer  (1976),  Thornton  and  Morris  (1978),  Downing  (1984), 
JafTe,  Sternberg,  and  Sallengcr  (1985),  Sternberg,  Shi,  and  Downing  (1989),  Hanes 
and  Huntley  (1986),  Beach  and  Sternberg  (1987),  and  Hanes  et  al.  (1988).  These 
observations  were  made  at  several  different  sites,  both  within  and  seaward  of  the 


941 


942 


D.  M.  HANES 


surf  zone,  and  using  several  different  instruments.  A  brief  discussion  of  some  of 
these  observations  within  the  context  of  intermittent  suspension  follows. 

In  the  region  just  seaward  of  the  wave  breaking  zone.  Hanes  and  Huntley  (1986) 
deployed  a  miniature  optical  backscatter  sensor  (MOBS)  at  Pt.  Sapin.  New 
Brunswick,  and  observed  that  both  the  suspension  of  sediment  and  the  cross-shore 
transport  of  suspended  sediment  were  intermittent.  The  suspended  sediment  was 
associated  with  incident  waves  and  was  enhanced  during  wave  groups.  Within  the 
surf  zone,  observations  of  intermittent  suspension  have  been  made  using  fast 
response  optical  instrumentation.  Using  light  transmission.  Brenninkmeyer  (1976) 
observed  infrequent  bursts  of  high  concentrations  ( >  10  g/1)  of  suspended  sediment 
in  the  surf  zone  and  frequent  bursts  of  extremely  high  concentrations  (up  to  500 
g/1)  of  suspended  sediment  in  the  swash  zone.  Downing  et  al.  (1985).  Jaffa. 
Sternberg,  and  Sallenger  (1985),  and  Beach  and  Sternberg  (1987)  have  used  optical 
backscatter  sensors  to  observe  intermittent  suspension  within  the  surf  zone,  which 
exhibited  time  scales  corresponding  to  both  incident  and  infragravity  wave  periods. 
Jaffe,  Sternberg,  and  Sallenger  (1985),  in  particular,  point  out  the  significance  of 
intermittent  suspension  to  cross-shore  sediment  transport  by  showing  the  signifi¬ 
cance  of  the  flux  coupling  between  fluid  velocity  and  sediment  suspension. 


2.  The  Stanhope  Lane  Experiment 

The  experiment  al  Stanhope  Lane.  Prince  Edward  Island.  Canada,  in  October.  1984 
was  a  component  of  a  larger  study  designed  to  e.xamine  some  of  the  nearshore 
processes  that  occur  during  storms  and  to  help  understand  how  barred  beaches 
respond  to  storms.  Gillie  (1985)  and  Daniel  (1986)  provide  detailed  descriptions  of 
the  overall  C2S2  project.  Descriptive  obser%'aiions  of  suspended  sediment  obtained 
with  an  acoustic  concentration  meter  (ACM)  have  been  resorted  in  Hanes  ei  al. 
(1988). 

The  acoustic  system  that  was  used  to  obtain  the  measurements  to  be  described 
was  constructed  for  the  NO.A.A  .Atlantic  Oceanographic  and  Meteorological  Labora¬ 
tories  in  Miami,  Florida,  and  was  first  described  by  Huff  and  Fiske  (1980).  The 
primary  advantages  of  the  .ACM  are  that  it  is  nonintrusive  and  with  suitable 
calibration  can  measure  concentration  profiles  with  high  resolution.  The  initial  tank 
calibrations  were  carried  out  by  Young  et  al.  (1982),  with  more  detailed  calibration 
for  the  C2S2  experiment  described  in  Hanes  et  al.  (1988).  The  data  to  be  presented 
nave  temporal  resolution  of  0.8  sec  and  spatial  resolution  of  approximately  10 
percent  of  the  distance  above  the  sea  bed  for  any  particular  measurement,  but  not 
les.s  than  0.5  cm.  The  return  due  to  the  sea  bed  biases  the  concentration  measure¬ 
ments  within  approximately  1.5  cm  of  the  sea  bed. 

The  gain  of  the  .ACM  was  chosen  to  avoid  saturation  at  high  concentrations.  .At 
the  gain  setting  used  for  this  experiment,  the  minimum  resolution  in  concentration 
measurements  was  approximately  70  mg/'l.  .Although  this  represents  a  high  sus¬ 
pended  sediment  concentration  in  deeper  water,  the  concentration  at  the  measure¬ 
ment  site  during  times  of  interest  was  on  the  order  of  1000  to  5000  mg/1.  The  gain 
was  therefore  chosen  to  enable  measurement  of  concentration  in  the  100  to  10.000 
mg/1  range. 

The  data  to  be  discussed  in  this  chapter  were  obtained  at  a  water  depth  of 
approximately  2  m.  near  the  crest  of  a  longshore  sand  bar,  designated  as  position 
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Fig.  1.  Beach  profile  at  experimental  site  showing  location  (number  2)  of  instrument  deployment 


number  2  in  Fig.  1.  The  sediment  consisted  of  well  sorted  sand  with  a  modal  sieve 
size  of  0.21  mm.  Instrumentation  consisted  of  an  acoustic  concentration  meter  that 
probed  the  lowest  meter  of  the  water  column  and  an  electromagnetic  flowmeter 
located  10  cm  above  the  sea  bed.  The  measurements  were  obtained  over  a  410-sec 
period  following  a  several  day  long  severe  storm.  The  tide  was  falling  and  sea  level 
was  near  its  mean  value.  The  peak  incident  wave  period  was  6,1  sec,  the  cross-shore 
rms  fluid  speed  as  measured  by  the  flowmeter  was  24.0  cm/sec,  and  the  rms  wave 
height  is  estimated  to  have  been  24  cm.  We  were  unable  to  obtain  accurate 
quantitative  bed  form  measurements  during  the  measurement  period,  but  visual 
observations  indicated  ripples  with  lengths  of  order  10  cm.  This  data  set  corre¬ 
sponds  to  set  D  in  Hanes  et  al.  (1988). 

3.  Structure  of  Events 

The  time  averaged  suspended  sediment  concentration  profile  is  given  in  Fig  2.  The 
mean  concentration  increases  sharply  in  the  7-cm  region  near  the  sea  bed.  This  is 
the  region  of  interest  for  the  remainder  of  this  paper.  The  depth  average  of  the 
instantaneous  suspended  sand  concentration  between  the  limits  of  2  and  7  cm  above 
the  sea  bed  will  be  referred  to  in  what  follows  as  the  near  bottom  concentration 
(NBC). 

A  lime  series  of  the  NBC  and  fluid  speed  U  is  shown  in  Fig.  3,  where  positive  U 
indicates  that  the  cross-shore  component  of  the  fluid  velocity  is  directed  offshore. 
Tne  suspended  sediment  concentration  is  clearly  unsteady  and  dominated  by  events. 
The  temporal  structure  and  spatial  siruaure  of  these  suspension  events  are  clarified 
by  examining  the  variations  in  concentration  on  short  time  and  space  scales.  We  will 
first  examine  a  particular  event,  followed  by  a  composite  description  of  all  the 
events. 

A  particular  event  that  occurred  approximately  in  the  middle  of  the  record  will 
be  described.  A  time  series  of  concentration  at  various  elevations  (all  less  than  7  cm) 
above  the  sea  bed  and  fluid  speed  is  shown  in  Fig.  4,  where  the  elapsed  time 
between  measurements  is  0.8  sec.  The  instantaneous  concentration  during  the  event 
decreases  with  elevation  above  the  sea  bed.  as  expected  from  the  mean  concenira- 
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Fig.  2.  The  vertical  profile  of  suspended  sediment  concentration  averaged  over  410  sec. 


tion  profile.  This  tendency  is  prevalent  in  all  events.  This  suspension  event  has  a 
time  scale  of  approximately  3  to  4  sec  and  a  vertical  spatial  scale  of  approximately 
5  cm.  For  comparison,  the  settling  velocity  of  the  sand  JVq  is  approximately 
2  cm/sec  and  U  during  the  event  averaged  approximately  -20  cm/sec. 

A  composite  description  of  the  intermittent  events  is  obtained  through  a  condi¬ 
tional  sampling  technique.  By  setting  a  threshold  equal  to  the  mean  plus  one 
standard  deviation  of  the  NBC  for  the  entire  record.  20  events  were  identified.  These 
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Fig.  3.  Time  series  of  depth  averaged  (between  2  and  7  cm  above  the  sea  bed)  suspended  sediment 
concentration  (NBQ  and  fluid  speed,  where  positive  U  indicates  the  cross-shore  component  of  the  fluid 
speed  is  directed  offshore. 
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Fig.  4.  Time  senes  of  suspended  sediment  concentration  at  several  elevations  and  fluid  speeds,  for  a 
particular  event  centered  near  pulse  269. 


Time,  sec 


Fig.  S.  Composite  average  for  20  suspension  events  of 
suspended  sediment  concentration  at  2  to  6  cm  elevation 
in  ^-cm  increments.  The  concentration  decreases  mono- 
tonically  with  elevation,  with  the  upper  curve  corre¬ 
sponding  to  2  cm  above  the  sea  bed  and  the  lowest  curve 
at  6  cm  above  the  sea  bed. 


20  events  and  the  corresponding  velocity  measurements  are  then  ensemble  averaged 
by  temporally  aligning  their  maximum  concentrations.  The  resulting  average  event 
is  shown  in  Fig.  5.  This  composite  event  lasts  for  about  3  to  4  sec  and  has  a  peak 
concentration  (at  elevation  2  cm)  of  about  4.5  g/1.  The  concentration  decreases 
sharply  with  distance  above  the  sea  bed.  with  the  concentration  variations  negligible 
above  about  6  cm.  .\iso  notable  is  the  rapid  increase  and  subsequent  rapid  decrease 
in  concentration  with  time. 


4.  Occurrence  of  Events 

The  cause  of  the  intermittenc)'  of  these  near  bed  suspension  events  is  one  of  the 
perplexing  issues  that  immediately  arises  from  these  observations.  There  are  two 
interesting  observations  of  relevance  to  this  question.  First,  the  occurrence  of  events 
is  correlated  with  periods  of  higher  velocity  associated  with  groups  of  larger  waves; 
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Fig.  6.  Cross  spectral  coherence  between  the  NBC  and  SKE  indicates  the  correlation  between  wave 
groups  and  intermittent  suspension  events  near  periods  of  30  sec.  Dashed  lines  correspond  to  U,  solid 
lines  to  SKE.  and  the  lowest  solid  line  of  the  upper  plot  is  the  autospectrum  of  NBC. 


second,  they  initiated  (in  16  out  of  20  occurrences)  during  the  shoreward  phase  of 
the  cross-shore  fluid  motion. 

From  Fig.  3  it  can  be  seen  that  the  suspension  events  are  usually  in  the  midst  of 
wave  groups  and  rarely  occur  coincident  with  the  first  wave  of  the  group.  The 
groupiness  in  the  velocity  record  is  quantified  by  the  smoothed  kinetic  energy 
(SKE),  denned  as  SKE(/)  =  where  is  the  peak  incident  wave 

period  and  U  is  the  velocity  time  series.  The  correlation  between  the  NBC  and  the 
SKE  is  0.4,  which  is  significant  with  above  a  99  percent  probability.  The  cross 
spectral  coherence  between  the  NBC  and  the  SKE  is  shown  in  Fig.  6.  Although  the 
short  duration  of  the  overall  measurement  period  relative  to  the  group  time  scale 
results  in  large  confidence  intervals,  the  peak  in  coherence  near  a  period  of  30  sec  is 
interpreted  as  indicative  of  the  existence  of  a  relationship  between  wave  groups  and 
suspension  events  (which  is  not  yet  well  described). 

A  particularly  striking  example  of  an  isolated  suspension  event  that  occurred  in 
the  middle  of  a  wave  group  (at  approximately  130  sec  in  Fig.  3)  is  shown  with  an 
enlarged  scale  in  Fig.  7.  Tnis  particular  event  is  of  very  short  duration  and  occurs 
during  the  shoreward  phase  of  motion  corresponding  to  the  largest  wave  in  the 
group.  It  is  likely  that  there  is  a  threshold  for  these  events  that  is  exceeded  by  either 
the  larger  waves  within  a  group  or  some  phenomenon  associated  with  the  sequence 
of  larger  waves. 

Over  the  entire  record  length,  the  correlation  between  U  and  the  rate  of  change  of 
NBC  was  0.26.  Although  significant,  this  correlation  only  explains  a  small  portion  of 
the  covariations  between  U  and  NBC.  The  relationship  between  the  entrainment 
of  sediment  from  the  sea  bed  and  the  cross-shore  fluid  velocity  is  illustrated  in  Fig.8. 
which  is  a  scatter  diagram  of  the  rate  of  change  of  NBC  verses  U,  for  times  when  the 


INTERMITTENT  SUSPENSION  OF  SAND 


947 


®?To  TTs  120  125  130  135  UO  US  150 

Time,  s 


Fig.  7.  Enlarged  section  of  Fig.  3  illustrating  an  isolated  suspension  event  that  occurred  in  the  midst  of 
a  group  of  larger  waves. 


NBC  is  increasing.  The  periods  of  rapidly  increasing  NBC  occur  more  frequently 
while  the  cross-shore  velocity  is  directed  onshore  than  when  the  velocity  is  directed 
offshore.  Because  the  peaks  in  onshore  velocity  were  only  slightly  greater  than  the 
peaks  in  offshore  velocity,  this  phenomenon  cannot  be  explained  in  terms  of  a 
simple  magnitude  threshold.  Comparing  periods  of  equivalent  magnitude  onshore 
and  offshore  velocity,  for  example  for  40  <  |  C/J  <  50  cm/sec,  the  NBC  during  the 
onshore  motion  is  a  factor  of  3.3  greater  than  during  the  offshore  motion.  Qualita¬ 
tively  similar  observations  were  also  made  by  Hanes  and  Huntley  (1986)  at  a 
different  field  site  using  different  instrumentation. 

Hallermeier  (1982),  Sleath  (1982),  and  Hanes  and  Huntley  (1986)  have  suggested 
that  fluid  acceleration  may  be  causally  related  to  the  concentration  of  suspended 
sand.  For  the  present  measurements,  the  correlation  between  water  acceleration  and 
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Fig.  8.  Scatter  plot  of  the  rate  of  change  of  the  NBC  versus  U 
for  times  when  NBC  was  increasing. 
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Fig.  9.  Scatter  plot  of  the  rate  of  change  of  NBC  versus  the  offshore  fluid  acceleration. 


NBC  was  0.084,  which  is  barely  significant  at  the  95  percent  level.  However,  the 
anticorrelation  between  acceleration  and  the  time  rate  of  change  of  NBC  is  about 
twice  as  high  (  —  0.164),  as  shown  in  Fig.  9.  This  supports  the  idea  that  water 
acceleration  is  significant  to  the  entrainment  of  sand  from  the  sea  bed,  but  the  low 
correlation  is  suggestive  that  the  connection  between  acceleration  and  sand  entrain¬ 
ment  is  moderated  by  other  processes  that  have  not  yet  been  quantified. 

5.  Cross-Shore  Sediment  Transport 

The  correlation  between  high  suspended  sediment  concentrations  and  shoreward 
fluid  motion  has  a  strong  influence  on  the  cross-shore  sediment  transport,  resulting 
in  a  net  shoreward  transport  of  sediment  during  the  measurement  period.  According 
to  linear  wave  theory,  there  is  little  vertical  structure  in  the  wave  induced  fluid 
motions.  The  suspended  sediment  flux  is  therefore  estimated  at  any  depth  as  the 
concentration  at  that  depth  times  the  cross-shore  fluid  velocity  measured  at  the 
current  meter  approximately  10  cm  above  the  sea  bed.  Huntley  and  Hanes  (1987) 
describe  a  method  to  correct  the  cross-shore  velocity  to  crudely  account  for  the 
bottom  wave  boundary  layer,  but  found  such  a  correction  had  little  influence  on  the 
trends  of  the  sediment  transport.  However,  the  results  to  be  presented  are  subject  to 
the  caveat  that  the  near-bed  water  velocity  structure  due  to  waves  over  a  rippled  bed 
have  not  been  well  documented  in  the  field;  there  may  be  structure  in  the  velocity 
field  in  the  bottom  10  cm  that  we  have  not  taken  into  account. 

The  NBC  and  estimated  cross-shore  suspended  sediment  flux  can  be  described 
using  the  composite  averaging  technique.  Figure  10  shows  the  results  for  the  20 
events,  where  the  four  curves  represent  the  composite  averages  of  the  NBC,  the 
NBC  plus  one  standard  deviation,  the  fluid  speed  U,  and  the  estimated  cross-shore 
suspended  sediment  flux.  The  event  initiates  as  the  cross-shore  velocity  transitions 
from  seaward  to  shoreward  and  the  peak  in  concentration  occurs  while  the  velocity 
is  directed  shoreward.  The  concentration  decreases  as  the  fluid  velocity  direction 
reverses  back  to  seaward  motion.  Although  the  average  cross-shore  velocity  is 
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Fig.  10.  Composite  average  of  20  events.  The  starred  carves 
represent  in  grains  per  liter  the  NBC  plus  one  standard 
deviation.  The  diamond  curve  is  U  in  meters  per  second  and 
the  triangle  curve  is  the  estimated  cross-shore  sediment  flux, 
in  grams  per  10  centimeters  per  centimeter  per  second. 


directed  seaward,  the  net  flux  during  the  event  is  directed  shoreward  because  of  the 
co-occurrence  of  the  peak  concentration  with  shoreward  velocity. 

Figure  11  shows  the  vertical  profile  of  cross-shore  sediment  flux  along  with  die 
profile  of  the  product  between  the  time  averaged  velocity  and  time  averaged 
concentration.  As  pointed  out  by  Jafie,  Sternberg,  and  Sallenger  (1985)  and  Hanes 
and  Vincent  (1987),  the  fluctuation  between  the  variations  in  concentration  and 
velocity  in  a  wave  dominated  environment  are  more  significant  to  cross-shore 
transport  than  the  product  of  the  mean  quantities.  Although  the  20  events  identified 
earlier  cumulatively  only  account  for  35  percent  of  the  record  length,  the  cross-shore 
sediment  transport  during  the  events  represents  82  percent  of  the  total  over  the 
entire  record,  .daringly,  the  20  points  corresponding  to  the  peak  concentrations  of 
the  20  events  account  for  59  percent  of  the  total  transport,  even  though  they  occur 
during  only  4  percent  of  the  total  measurement  period. 

Huntley  and  Hanes  (1987)  used  standard  time  series  analysis  techniques  to 
express  the  cross-shore  flux  in  terms  of  the  cospectrum  of  U  and  concentration.  It  is 
of  interest  to  apply  this  methodology  to  these  data,  even  though  the  transport  is 
dominated  by  events.  The  cospectrum  between  U  and  NBC  is  shown  in  Fig.  12.  Tne 
NBC  spectrum  is  fairly  broad,  but  the  cospectrum  exhibits  a  strong  negative  peak 
near  the  incident  wave  period,  which  corresponds  to  a  shoreward  transport  of 
suspended  sand. 
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Fig.  11.  Veriical  profile  of  Uie  time  averaged  cross-short 
suspended  sediment  flux  isquaresi  and  the  product  of  the  time 
averages  of  U  and  C  (iriangjesV  The  star  represents  the  depth 
and  time  averaged  flux. 
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Fig.  12.  The  cospectrum  of  U  and  NBC,  spectrum  of  U  (dashed)  and  spectrum  of  NBC  (solid).  A 
negative  cospectral  value  indicates  shoreward  transport. 


6.  Discussion 

There  are  multiple  processes  consistent  with  these  observations  that  are  suggestive 
of  rapid  vertical  entrainment  of  sediment  from  the  sea  bed.  The  precise  form  of  this 
entrainment  is  not  yet  known.  It  may  be  a  manifestation  of  the  two-dimensional 
ejection  and  advection  of  sediment  laden  vortices  due  to  the  rippled  bed,  as  seen  for 
example  in  the  laboratory  by  Sleath  (1982).  .Alternatively,  the  event  may  be  nearly 
one-dimensional,  consisting  of  the  vertical  motion  of  a  layer  of  sand,  as  modeled  for 
example  by  Fredsoe,  Andersen,  and  Silberg  (1985).  The  scales  of  the  events  support 
the  vertical  entrainment  scenario.  If  the  events  are  considered  to  be  the  result  of  the 
horizontal  advection  of  a  vortex  cloud  of  suspended  sediment  through  the  measure¬ 
ment  area,  then  the  duration  of  the  event  would  be  less  than  those  observed.  For 
example,  using  the  measured  vertical  length  scale  of  5  cm  as  the  length  scale  (I.)  of 
a  cloud  and  the  measured  average  horizontal  fluid  velocity  {U)  during  the  events  of 
6  cm/ sec.  then  the  time  scale  of  the  event  would  be  order  L/t/  =  0.83  sec,  which  is 
about  a  factor  of  5  too  low.  It  is  of  interest  to  note,  however,  that  there  were  events 
with  short  time  scales  such  as  this  on  other  days  (Pescio  and  Hanes,  1988).  If 
alternatively,  sand  were  entrained  from  the  sea  bed  and  then  settled  back  to  the  sea 
bed  at  a  velocity  W'q  of  2  cm/sec,  the  time  scale  would  be  approximately  1L/Wq  = 
5  sec,  which  is  consistent  with  the  observations. 

The  intermittency  of  the  events  observed  at  this  site  is  interesting  and  probably 
results  from  wave  conditions  that  are  near  the  threshold  for  sediment  entrainment 
The  sequence  of  waves  in  a  group  appears  to  enhance  suspension.  This  may  be  due 
to  a  hydrodynamic  phenomena,  which  indicates  that  the  vonicity  induced  near  the 
sea  bed  during  a  wave  half-cycle  may  not  be  fully  cancelled  by  the  oppositely 
directed  vonicity  induced  in  the  successive  half-cycle.  Alternatively,  these  observa- 


INTERMITTENT  SUSPENSION  OF  SAND 


951 


lions  could  result  from  the  internal  response  of  a  porous  bed  to  forcing  by  wave 
groups. 


7.  Concluding  Remarks 

The  classical  paradigm  of  suspended  sediment  in  unidirectional  flow  in  which  the 
tendency  of  the  sediment  to  sink  is  balanced  by  a  vertically  upward  diffusive  flux 
caused  by  fluid  turbulence  does  not  apply  in  wave  dominated  environments.  This  is 
not  surprising  because  many  investigators  have  described  the  organized  process  by 
which  waves  interact  with  a  rippled  bottom  to  form  and  eject  sediment  laden  eddies. 
Even  the  casual  swimmer  can  frequently  notice  the  formation  of  either  cloudy  bursts 
of  sediment,  or  a  rising  and  falling  sheet  of  sediment,  depending  on  the  wave  and 
sea  bed  conditions. 

Many  new  questions  follow-  from  the  observation  that  suspension  initiates  during 
the  shoreward  phase  of  the  fluid  motion  and  that  the  net  sediment  transport 
resulting  from  this  process  is  directed  shoreward.  For  example,  which  features  of  the 
forcing  are  responsible  for  this  asymmetry,  is  the  process  sometimes  reversed,  and 
does  it  occur  across  the  inner  continental  shelf?  Further  observations  of  suspended 
sediment  concentration,  sea  bed  microtopography,  and,  if  possible,  fluid  and  sedi¬ 
ment  velocity  at  high  resolution  near  the  sea  bed  and  within  the  sea  bed  would  help 
to  begin  to  answer  some  of  these  questions. 
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ABSTR.ACT 

Ludwig.  K..A.  and  Hanes.  D.M.,  1 990.  A  laboratory  evaluation  of  optical  backscatterance  suspended  solids  sensors  exposed 

to  sand-mud  mixtures.  Mar.  Geol.,  94:  173-179. 

Optical  suspended  solids  sensors  were  studied  in  the  laboratory-  in  order  to  obtain  calibration  factors  and  evaluate  the 
general  behavior  of  the  instrument  when  exposed  to  suspensions  of  mud  and  sand.  The  sensors  were  found  to  be  capable 
of  measunng  the  time-averaged  concentration  of  either  suspended  sand  or  suspended  mud.  but  the  relative  sensitivity  of 
the  sensor  to  these  two  different  types  of  sediment  varied  by  approximately  an  order  of  magnitude.  In  mixtures  with  a 
constant  or  known  mud  concentration,  the  instrument  response  to  the  mud  may  be  treated  as  an  offset,  and  gains  appro¬ 
priate  to  sand  may  be  used  to  estimate  the  suspended  sand  concentrations.  Semi-quantitative  observations  of  the  effects 
of  flow,  direction  relative  to  sensor  orientation  indicate  that  significant  concern  is  warranted. 


Introduction 

Optical  backscatterance  suspended  solids 
sensors  have  been  used  in  the  field  in  order  to 
study  the  dynamics  of  suspended  sediment  in 
coastal  regions  (Downing  et  al.,  1981;  Stern¬ 
berg  et  aJ..  1984;  Hanes  and  Huntley,  1986). 
These  sensors,  manufactured  by  D&.A  Instru¬ 
ments  and  Engineering,  emit  infrared  (IR) 
light,  and  measure  the  quantity  of  IR  backscat- 
tered  from  a  small  (order  of  10  cmM  concen¬ 
tration-dependent  volume  immediately  adja¬ 
cent  to  the  sensor  face.  Laboratoiy  calibration 
allows  conversion  of  this  backscattered  light 
intensity  into  suspended  sediment  concentra¬ 
tion  (Downing  and  Beach,  1989). 
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This  paper  describes  laboratory  tests  de¬ 
signed  to  calibrate  the  sensitivity  of  the  instru¬ 
ments  to  a  variety  of  suspended  sediments,  and 
also  to  investigate  the  general  behavior  of  the 
instruments.  One  of  the  important  behavioral 
characteristics  investigated  in  these  experi¬ 
ments  is  the  response  of  the  instruments  to 
sand  in  the  presence  of  suspended  mud.  These 
tests  were  prompted  by  field  investigations  of 
the  erosion  of  Louisiana  barrier  island 
shorefaces. 

Sediments 

Most  of  the  sediments  used  in  this  study  were 
collected  at  a  field  site  near  Fourchon.  Louisi¬ 
ana.  located  southwest  of  New  Orleans  near 
Grand  Isle.  Both  exposed  beach  sands  and  off¬ 
shore  bottom  samples  were  collected  at  this  lo¬ 
cation.  The  bottom  samples  were  collected  by 
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divers  from  approximately  2  km  offshore,  at 
an  instrument  deployment  site  in  approxi¬ 
mately  7  m  of  water.  Examination  of  these  bot¬ 
tom  samples  revealed  silt,  clay,  and  file  sand 
with  occasional  shell  fragments. 

Following  removal  of  the  coarse  shell  frag¬ 
ments  by  seiving,  size  analysis  indicated  that 
more  than  75%  of  the  sediment  was  less  than 
20  um  in  size.  .According  to  the  “Glossary  of 
Geology."  mud.  in  the  context  of  marine  geol¬ 
ogy.  is  “a  sticky,  fine-grained,  marine  detrital 
sediment,  either  pelagic  or  terrigenous”  (Bates 
and  Jackson,  1980).  As  this  definition  aptly 
describes  the  Fourchon  offshore  bottom  sam¬ 
ples.  throughout  the  remainder  of  this  paper 
these  samples  shall  be  referred  to  as  mud. 

Fourchon  sand  samples  were  collected  from 
the  beachface  on  the  exposed  barrier  island. 
The  sand  was  rinsed,  dried,  and  sieved  through 
the  following  sieve  stack  (jim ):  500,  355. 250, 
180.  125.  90.  63.  and  a  catch  pan.  Hereafter, 
the  size  of  the  sediment  will  reflect  the  sieve 
size  in  which  the  sediment  was  found. 

Sand  samples  taken  from  two  Florida  loca¬ 
tions  were  also  used  in  some  of  the  tests.  The 
first  of  these  locations  is  a  beach  approxi¬ 
mately  I  mi  north  of  Port  Canaveral,  near  the 
cit\  of  Cape  Canaveral,  Florida.  The  second 
Florida  location,  chosen  for  its  proximity  to  the 
laboratoiy-.  is  a  beach  on  the  Virginia  Key  side 
of  Bear  Cut  Inlet,  located  slightly  east  of 
Miami.  Florida.  Figure  1  shows  the  grain-size 
distributions  for  these  sand  sampies.  The  mean 
and  modal  seive  diameters  (um)  are  respec¬ 
tively.  1 9 1  and  1 87  for  Fourchon.  289  and  237 
for  Bear  Cut.  and  247  and  172  for  Cape 
Canaveral. 

Calibration  facility 

The  primary  test  apparatus  for  this  study 
consists  of  a  recirculation  calibration  tank  (Fig. 
2  )  containing  46  1  of  water  and  sediment  cir¬ 
culated  at  an  average  speed  of  2.2  cm/s  in  the 
test  area.  The  water  and  sediment  flow  down¬ 
ward  through  the  test  section  and  are  recircu- 


GRAIN  SIZE,  microns 
Fig.  I .  Beach  sand  size  distnbutions. 


Fig.  2.  Scheme  of  recirculation  calibration  facility,  includ¬ 
ing  water-sampling  apparatus. 

lated  with  a  centrifugal  pump.  The  discharge 
passes  through  a  diffuser  at  the  top  of  the  tank 
in  order  to  thoroughly  mix  the  sediment.  The 
test  section  of  the  tank  consists  of  a  vertical 
clear  acr>'lic  tube  with  ports  through  which 
water  samples  may  be  drawn.  These  ports  are 
located  30.  60.  and  90  cm  below  the  full  mark 
of  the  tank.  The  hydraulic  head  difference  be- 
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tween  the  top  of  the  tank  and  the  discharge  end 
of  the  clear  nylon  tube  drives  the  water-sam¬ 
pling  process. 

The  tank  was  designed  to  allow  for  the  crea¬ 
tion  of  a  uniform  and  known  concentration  of 
suspended  sediment.  No  attempt  was  made  to 
simulate  any  panicular  hydraulic  conditions 
which  might  be  found  in  a  field  situation.  Ef¬ 
fects  such  as  the  dependence  of  calibration 
coefficients  on  Reynold’s  number  will  not  be 
addressed,  although  concerns  related  to  flow 
direction  will  be  presented  later  in  this  paper. 

The  abiiit}'  to  create  suspensions  of  known 
concentration  was  verified  through  physical 
sampling  of  the  water-sediment  mixture. 
Samples  of  approximately  2  1  were  drawn 
through  the  \  in  tube  and  passed  through  a  63 
,um  sieve  w  hich  trapped  the  sediment.  Prior  to 
further  discussion  of  these  samples,  the  follow¬ 
ing  definitions  will  be  introduced: 

Input  concentration:  The  dry  weight  of  sedi¬ 
ment  placed  in  the  tank  (g)  divided  by  the  vol¬ 
ume  of  water  placed  in  the  tank  (1). 

Sampled  concentration:  The  dry  weight  of 
sediment  sampled  (g)  divided  by  the  volume 
of  water  sampled  (1). 

Actual  concentration:  The  dry  weight  of  sed¬ 
iment  suspended  in  the  test  section  (g)  di¬ 
vided  by  the  volume  of  water  (1). 

P:  The  ratio  of  the  sampled  concentration  to 
the  input  concentration. 

In  a  uniformly  mixed  suspension  with  no 
sediment  settling  onto  surfaces,  the  input  con¬ 
centration  is  equal  to  the  actual  concentration. 
This  is  indeed  the  case,  as  will  be  shown  below. 

When  water  samples  are  draw  through  a  tube 
oriented  perpendicular  to  the  mean  How  direc¬ 
tion.  the  resulting  sampled  concentration  is  ex¬ 
pected  to  be  slightly  lower  than  the  actual  con¬ 
centration.  In  a  similar  sampling  apparatus, 
Bosman  ( !  9SZ  i  found  that  when  the  sampling 
velocity  e.xceeds  the  mean  test  section  flow,  the 
ratio  of  sampled  concentration  to  actual  con¬ 
centration  approaches  a  value  of  0.73.  Water 
sample  velocities  obtained  in  the  present  tests 
varied  with  the  elevation  of  the  discharge  end 


of  the  sampling  tube.  At  a  level  of  30  cm  below 
the  initial  tank  level,  a  sample  velocity  of  85 
cm/s  is  obtained.  At  90  cm  below  the  initial 
tank  level,  the  sample  velocity  increases  to  1 60 
cm/s.  In  all  cases  the  velocities  greatly  exceed 
the  cross-sectionally  averaged  flow  velocity  of 
2.2  cm/s. 

Tests  were  conducted  to  assess  the  possible 
dependence  of  the  sampled  concentration  on 
the  venical  location  of  the  sample  intake,  the 
cross-sectional  location  of  the  sample  intake 
(1),  the  input  concentration  of  sediment,  and 
the  flow  rate.  The  results  of  all  these  tests  ap¬ 
pear  in  Tables  1-3.  None  of  these  variables  has 
a  significant  and  consistent  effect  on  the  P 
coefficient.  The  fact  that  P  is  about  0.85,  rather 
then  0.73  as  found  by  Bosman,  is  probably  due 
to  the  flow  intake  geometry.  However,  the 
value  of  P  is  not  significant  to  the  reminder  of 
the  results;  these  tests  were  simply  performed 
to  verify  that  the  tank  was  well  mixed.  There- 

TABLE  ! 


ji  versus  venical  location  and  sample  flow  rate 


Sample  head 

30  cm 

60  cm 

90  cm 

Pon  30 

0.S3 

0.95 

0.80 

Port  60 

0.93 

0.83 

0.88 

Pon  90 

0.88 

0.85 

0.85 

Input  concentration  =  1 .0  g/1. 

Port  number  indicated  distance  ( cm )  below  full  level. 
Intake  at  center  of  tank. 


TABLE ; 

P  versus  input  concentration 


Cone,  (g/1)  p 


0.05 

0.55 

0.10 

0.80 

0.50 

o.s: 

1.0 

0.85 

5.0 

0.-9 

10.0 

0.78 

Port  90. 

Head  90  cm. 
Intake  /=9..'  cm. 
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TABLE  3 


/I  versus  / 

/  (cm ) 

a 

1.27 

0.S9 

9.50 

0.S5 

9.50 

0.S3 

14.60 

0.S5 

Port  90. 

Head  90  cm. 

Input  Cone.  =  1.0  g/1. 


fore,  the  input  concentration  may  be  assumed 
to  be  equal  to  the  actual  concentration.  This  is 
believed  to  be  a  valid  assumption  for  all  the 
tests,  independent  of  sediment  type  or  size,  and 
is  a  desirable  feature  of  the  recirculation  tank. 

Because  the  fluid  in  front  of  the  sensor  face 
is  in  constant  motion,  it  is  reasonable  to  as¬ 
sume  that  there  will  be  a  variation  in  concen¬ 
tration  with  respect  to  time.  Physical  samples 
indicate  that  when  this  variation  is  averaged 
over  appro.ximately  1  min.  the  resulting  sam¬ 
pled  concentration  becomes  steady,  and  rep¬ 
resentative  of  the  input  concentration.  Simi¬ 
larly,  time  averaging  of  the  OBS  signal  also 
yielded  a  steady  mean  voltage,  which  may  then 
be  equated  with  the  input  concentration. 

Results  and  discussion 

In  order  to  test  the  response  of  OBS  sensors 
to  suspended  beach  sand,  measurements  were 
taken  using  the  three  different  sands  at  several 
concentrations  ranging  from  clear  water  to  10.0 
g/1.  In  clear  tap  water  the  signal  fluctuations 
w'ere  approximately  equal  to  the  resolution  of 
the  analog-To-digital  convener,  about  5  mV. 
.A.lthough  the  gains  and  offsets  vary  between 
sensors,  each  sensor  is  quite  linear  over  the 
range  of  concentrations  considered.  Each  sen¬ 
sor  has  an  adjustable  gain,  so  the  difference  be¬ 
tween  sensor  gains  is  superfluous  to  this  dis¬ 
cussion.  The  correlation  coefficients  typically 
exceed  0.99.  .An  example  of  the  results  of  these 
tests  for  Fourchon  beach  sand  appears  in 


graphical  form  in  Fig.  3.  The  complete  list  of 
gains  is  given  in  Table  4. 

An  important  consideration  in  the  deploy¬ 
ment  of  these  sensors  at  a  field  site  which  con¬ 
tains  sands  of  varying  sizes  is  the  dependence 
of  the  response  on  the  sediment  size.  The  re¬ 
sponse  to  different  grain  sizes  was  e.xamined 
by  seiving  the  Fourchon  beach  sand  into  size 
classes  of  63.  90,  125,  180  and  250  urn  and  a 
mixture  which  approximated  the  natural  dis¬ 
tribution.  The  mean  grain  size  of  the  mixture 
was  183  um. 

Figure  4  shows  the  gain  versus  mean  size  of 
each  size  fraction  for  sensors  1  and  4.  For  each 
sensor,  gains  for  sands  ranging  in  size  from  63 
to  355  /tm  differ  from  the  gain  of  the  mixture 
by  a  maximum  of  about  17%.  The  variation 
with  size  is  not  monotonic,  presumably  be- 


Fig.  3.  Response  to  Fourchon  beach  sand. 


T.ABLE  ■i 


Gains  for  the  \  anous  sediments 


Sieved 

mud 

Fourchon 

sand 

Canaveral 

sand 

Bear  Cut 
sand 

Sensor  ! 

6s6.! 

"6.5 

93.5 

98.5 

Sensor  2 

115.1 

12.8 

15.5 

21.8 

Sensor  3 

^22 

67.6 

84.8 

84.7 

Sensor  4 

484.6 

50.4 

69.4 

71.9 

Sensor  5 

594,5 

50.5 

69.3 

78.4 
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Fig.  4.  Gain  of  sensors  I  and  4  as  a  function  of  sand  size. 

cause  of  the  differing  mineralogies  of  the  var¬ 
ious  size  fractions:  the  63  pim  size  group  ap¬ 
pears  nearly  black  in  color,  suggesting  high 
concentrations  of  magnetite  and  hornblende, 
while  the  tan  coloration  of  the  groups  over  125 
,um  suggests  compositions  of  largely  quartz  and 
feldspar. 

Two  series  of  tests  were  conducted  using 
suspended  mud.  The  first  series  determined 
calibration  gains  and  offsets  for  the  range  of 
concentrations  anticipated  in  the  field.  The 
concentrations  of  mud  used  in  these  tests  were 
0.0.  0.45,  0.9  and  1.8  g/1.  Results  for  each  sen¬ 
sor  appear  in  graphical  form  in  Fig.  5  and  in 
tabular  form  in  the  “sieved  mud"  column  of 
Table  4.  .4  second  series  of  tests  investigated 
the  behavior  of  these  instruments  over  their 
entire  measurement  range.  These  concentra¬ 
tions  were  staned  at  approximately  2  g/1  and 
increased  in  2  g/1  steps  until  sensor  one.  the 
most  sensitive  sensor,  registered  nearly  the 
same  voltage  for  two  successive  steps.  The  fi¬ 
nal  concentration  was  approximately  10  g/1. 
Figure  6  contains  the  results  of  this  series  of 
runs  for  all  five  sensors.  The  instrument  re- 


Fig.  5.  Response  of  all  five  sensors  to  moderate  concentra¬ 
tions  of  Fourchon  mud. 


Fig.  6.  Response  of  all  five  sensors  to  high  concentrations 
of  Fourchon  mud. 

sponse  to  mud  is  linear  up  to  concentrations  of 
approximately  3  or  4  g/1.  These  gains  are  nearly 
ten  times  as  great  as  those  determined  previ¬ 
ously  for  sand.  At  higher  concentrations,  the 
sensitivity  decreases,  although  not  markedly. 
Finally,  at  high  concentrations  which  depend 
on  the  resistor  settings  for  the  particular  sen¬ 
sor,  the  sensor  will  saturate  and  become  insen¬ 
sitive  to  funher  increases  in  concentration. 

In  order  to  recreate  a  condition  similar  to 
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that  expected  at  the  Louisiana  field  site,  var¬ 
ious  sand  concentrations  are  studied  in  two 
different  mud  backgrounds.  The  concern  is  the  | 
ability  of  the  sensors  to  “see”  the  sand  through  * 

a  background  mud  signal.  The  background  § 

mud  concentrations  have  values  of  0.0.  0.9  and  = 

o 

1.8  g/1.  Tests  are  run  on  these  backgrounds  = 
alone  and  with  these  backgrounds  in  combi-  | 
nation  with  concentrations  of  Fourchon  beach 
sand.  The  sand  concentrations  are  0.5.  1.0.  2.0 
and  5.0  g/1. 

The  results  of  these  runs  appear  in  Fig.  7  and 
in  Table  5.  The  “change  %”  column  of  Table  5 
contains  the  percentage  change  of  the  preceed- 
ing  gain  as  compared  to  the  gain  for  zero  mud 


SAND  CONCENTRATION,  g  i 

Fig,  Response  of  sensor  1  lo  sand  in  the  presence  of 
various  mud  bac.kgrounds. 


TABLE  5 


Gains  for  Fourchon  sand  in  various  mud  backgrounds 


Mud  =  0  g/1 

Mud  = 

0.9  g/1 

Mud  = 

l.Sg/l 

Gain 

Gain 

Change  “o 

Gam 

Change  °6 

Senior  1 

"6.55 

'6.50 

-O.i 

~6.20 

-0.5 

Sensor  2 

i;.s3 

1 1.62 

-9,4 

10.84 

-15.5 

Sensor  3 

67.62 

59.20 

-ll.S 

65.04 

-3.8 

Sensor  4 

50.44 

46.14 

-S.5 

44. '7 

-11.2 

Sensor  5 

50.52 

50.40 

-0.2 

50. '4 

-0,4 

1.0 

0.9 


o.sh  o 

o 

,  j, _ I _ I _ I _ I _ 

“■^0  AS  90  135  180  225 

RELATIVE  FLOW  DIRECTION,  degrees 

Fig.  8.  Sensor  1  output  as  a  function  of  sensor  orientation 
relative  to  the  mean  flow  direction.  .Angie  of  zero  corre¬ 
sponds  to  flow  directly  toward  the  sensing  face  of  the 
instrument. 

background.  The  gains  of  sensors  1  and  5  re¬ 
main  constant  to  within  1  %  regardless  of  the 
mud  background.  Sensors  2  and  4,  however, 
show  drops  in  gain  of  the  order  of  1 0  %  for  both 
mud  backgrounds  as  compared  to  sand-only 
concentrations.  Stranger  still  is  the  behavior  of 
sensor  3.  The  gain  experienced  by  this  sensor 
drops  11.8%  for  the  first  mud  background  and 
recovers  to  within  3.8  %  for  the  higher  mud 
background.  We  have  no  explanation  for  the 
difference  in  behavior  between  the  sensors.  Al¬ 
though  these  results  are  not  perfectly  consis¬ 
tent,  they  do  show  that  mud  backgrounds  of 
this  magnitude  reduce  the  gain  by  factors  in  the 
range  of  10  %  or  less. 

A  series  of  e.xperiments  examined  the  sensor 
response  as  a  function  of  the  orientation  of  the 
sensor  relative  to  the  direction  of  the  flow.  For 
these  experiments.  5.0  g/1  of  Bear  Cut  sand  was 
input  and  sensor  1  was  positioned  25  cm  below 
the  water  surface.  The  diffusion  nozzle  was  re¬ 
moved  from  the  discharge  hose  in  an  attempt 
to  obtain  a  high-energy  unidirectional  flow 
field.  No  sand  was  observed  accumulating  on 
any  part  of  the  sensor,  regardless  of  the  orien¬ 
tation  of  the  sensor.  Figure  8  indicates  that  the 
measured  concentration  varies  with  sensor 
orientation.  The  magnitude  of  the  variations 
reach  as  much  as  50%.  The  results  are  believed 
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to  result  from  the  flow  pattern  around  the 
sensor. 

Conclusions 

The  following  conclusions  may  be  drawn. 
The  sensors  possess  a  very  linear  response  for 
concentrations  of  sand  between  0.0  and  at  least 
10  g/1.  and  mud  between  0.0  and  2.0  g/1.  The 
sensor  sensitivity  to  suspended  mud  is  signifi¬ 
cantly  higher  than  that  due  to  sand.  .At  high 
mud  concentrations,  a  non-linear  calibration 
curve  is  required.  For  mixtures  of  mud  and 
sand,  the  sensor  responds  linearly  to  the  con¬ 
centration  of  each  component,  with  the  gains 
for  such  mixtures  being  comparable  to  those 
which  would  apply  to  each  component  alone. 

The  size  of  the  sand  in  suspension  does  not 
greatly  afTect  the  gain  of  the  sensors.  This  is  true 
for  the  grain-size  distribution  and  related  com¬ 
position  of  the  sand  found  at  the  Fourchon 
beach  site.  .A  more  consistent  dependence  of 
gain  on  size  might  be  expected  in  a  more  mi- 
neralogically  uniform  sand  sample. 

Use  of  these  instruments  for  the  measure¬ 
ment  of  suspended  sand  in  areas  which  con¬ 
currently  experience  suspended  mud  is  not 
recommended,  for  two  reasons.  First,  in  order 
to  avoid  saturating  the  sensors  with  mud  dur¬ 
ing  high-energy  events  ( these  are  typically  the 
events  of  most  interest),  the  instrument  gain 
must  be  set  so  low  that  high  resolution  of  sand 
suspension  is  prevented.  Second,  dunng  times 
w'hen  the  mud  concentration  vanes,  it  is  un¬ 
clear  how  a  distinction  should  oe  drawn  be¬ 
tween  the  ponion  of  the  return  due  to  sus¬ 
pended  sand  and  the  portion  of  the  return  due 
to  suspended  mud. 

•Although  the  flow  direction  dependence  tests 
do  not  attempt  to  duplicate  conditions  ex¬ 


pected  in  the  field,  they  do  indicate  that  a  sig¬ 
nificantly  strong  dependence  probably  exists. 
This  is  especially  true  in  the  wave-dominated 
environments  of  the  nearshore,  where  the  flow 
direction  is  frequently  changing.  For  deploy¬ 
ment  in  wave-dominated  environments  it  is 
recommended  that  the  sensors  be  oriented  such 
that  the  wave-induced  flow  is  parallel  to  the 
sensor  face.  Funher  testing,  preferably  in  a 
larger  tank  with  a  more  controlled  flow  struc¬ 
ture.  is  highly  recommended. 
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ABSTRACT 

Vincent.  C.E..  Hanes.  D.M.  and  Bowen.  A.J..  1991.  .Acoustic  measurements  of  suspended  sand  on  the  shoreface  and  the 
control  of  concentration  by  bed  roughness.  .Mar.  Geol..  96:  1-18. 

Acoustic  backscatter  measurements  of  the  concentrations  of  sand  in  s_.pe..sion  on  ihe  shoreface.  seawards  of  the  breaking 
zone,  during  a  mild  storm  event  show  that  sand  concentrations  incre..je  initially  but  then  rapidly  decrease  as  the  wave  energy 
increases:  it  is  suggested  that  the  bed  roughness  is  a  maio;  control  on  suspended  sand  concentration  and  that  the  decrease 
m  concentration  is  due  to  decreasing  rippie  steepness  after  the  break-off  bed  shear  stress  is  e.vtceded.  No  direct  measurements 
of  bedforms  were  available  and  the  combined  wave-current  interaction  model  of  Grant  and  Madsen  (J.  Geophys.  Res..  Vols. 
S-J  ( 1979)  and  87  (1982))  with  extensions  to  movable  bed  roughness  (ripple  dimensions  from  Carstens  et  al's.  1969  laboratory- 
data  (U.S.  Army  Corps  Eng.  CERC  Tech.  Mem.  28)  and  to  sediment-induced  seif-stratification  (Glenn  and  Grant.  J. 
Geophys.  Res..  Vol.  92,  1983)  was  used  to  predict  the  current  and  suspended  sand  concentration  profile.  Matching  the 
measured  and  modeled  concentrations  at  2  cm  above  the  bed.  the  values  of  the  resuspension  coefficient  Vj,  were  found  to 
decrease  as  excess  skin  friction  increased,  in  a  manner  similar  to  that  suggested  by  Drake  and  Cacchione  (Cent.  Shelf  Res.. 
Voi.  9.  1989)  but  were  an  order  of  magnitude  larger.  Using  the  much  smaller  rippie  dimensions  of  Nielsen  (J.  Geophys.  Res.. 
\  oi.  86.  1981)  changed  the  mode!  results  very  little.  Significant  differences  were  found  between  the  time-averaged  suspended 
sand  profiles  and  those  predicted  by  the  model. 

Sand  transport  fluxes  have  also  been  computed  using  the  instantaneous  products  of  measured  suspended  concentrations 
and  currents  (the  current  being  the  sum  of  steady  and  wave-induced  components).  The  wave-induced  fluxes  .show  considerable 
vanabiiity  from  run  to  run  but  the  general  pattern  is  of  a  transport  profile  having  a  shoreward  maximum  at  5-10  cm  above 
the  bed  and  offshore  transport  in  the  few  centimetres  closest  to  the  bed:  little  net  transport  occurred  above  15  cm.  These 
fluxes  show  a  weak  dependence  on  the  wave  energy,  becoming  more  shoreware  and  the  height  of  the  maximum  shoreward 
transport  decreasing  as  the  wave  energy  increased,  again  consistent  with  the  influence  of  the  rippie  steepness. 


1  Introduction 

Most  of  us  are  familiar  with  the  temporal 
changes  which  occur  on  beaches.  Storms  and  other 
periods  of  high  wave  activity  tend  to  move  sand 
offshore;  beach  levels  drop  and  the  land  behind 
the  beach  becomes  more  vulnerable  to  wave  attack, 
to  erosion  or  to  flooding.  Dunng  fair  weather, 
often  referred  to  as  "summer  wave  conditions",  a 
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gradual  accretion  occurs  and  beach  levels  rise. 
These  shorter  term  changes  are  often  accompanied 
by  more  gradual  changes  in  beach  volume  or 
position  which  may  be  due  to  alterations  in  sand 
supply  to  the  beach  (natural  or  man-induced)  or 
to  changes  in  the  wave  climate. 

Such  qualitative  descriptions  of  how  beaches 
alter  mask  the  problems  which  e.xist  in  quantifying 
the  processes  which  are  responsible  for  such 
changes.  The  alongshore  transport  of  sand  driven 
by  waves  breaking  on  the  shoreface  at  an  angle 
can  result  in  changes  in  beach  levels  if  there  is  a 
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flux  divergence  along  the  beach.  However,  most 
of  the  rapid  changes  in  beach  volumes  are  due  to 
the  cross-shore  transport  of  sand;  sand  is  trans¬ 
ported  from  the  upper  beach  face  out  to  bars 
beyond  the  breaker  zone  and  beach  levels  may  be 
reduced  by  many  tens  of  centimetres  in  the  course 
of  a  single  storm. 

The  data  described  here  were  collected  as  part 
of  the  C-CO.AST  field  experiment  designed  to 
investigate  the  hydrodynamic  and  sediment  trans¬ 
port  processes  on  beaches.  The  Canadian  Coastal 
Sediment  Transport  Program  (C-CO.AST)  is  a 
combined  project  of  the  Universities  of  Toronto. 
Dalhousie  and  Memorial  and  is  sponsored  by  the 
National  Science  and  Engineering  Research  Coun¬ 
cil  of  Canada.  The  results  presented  here  are  from 
a  single  acoustic  concentration  meter  (ACM)  and 
two  electromagnetic  current  meters  mounted  close 
to  the  sea  bed  100  m  from  the  mean  water  mark 
during  a  fieid  experiment  on  Queensland  beach. 
Nova  Scotia,  during  October  1987.  Results  from 
the  other  instrumentation  recording  concurrently 
will  be  presented  elsewhere. 

2  Queensland  beach  and  the  experimental  design 

Queensland  beach  is  a  small  sand\  pocket  beach 
some  300  m  wide,  enclosed  by  rocky  headlands. 
To  the  south  it  opens  into  a  larger  ba\  which  itself 
opens  to  the  .Atlantic  Ocean  (Fig.l).  .Atlantic  swell 


n;"  OuL-jnsiand  nejch.  Nciva  Scotia, 


and  storm  waves  penetrate  through  the  bay  open¬ 
ing,  arriving  as  long-crested  waves  almost  norma! 
to  the  beach.  The  beach  profile  is  a  simple  one 
with  an  upper  shore  face  slope  of  5'  at  mean  water 
level  (MWL).  decreasing  to  1.7°  100  m.  offshore 
(Fig.2).  .A  surncial  sand  sample  collected  by  divers 
100  m  offshore  close  to  the  .ACM  had  a  modal 
size  of  between  2.5  and  3.0®  when  dry  sieved 
through  half-phi  sieves  (Fig. 3). 

Instrumentation,  consisting  of  electromagnetic 
current  meters,  optical  backscatter  probes  (OBS). 
acoustic  concentration  meters  (ACM)  and  pressure 
sensors,  was  mounted  on  bottom  frames  60.  70. 
90  and  100  m  from  MWL  (Fig.2).  All  were  shore 
connected  and.  with  the  exception  of  the  ACMs. 
the  data  were  recorded  on  the  Dalhousie  University 
UDATS  system  (Hazen  et  ah.  1988).  Because  of 
their  high  data  return  rate  the  .ACMs  recorded 
directly  on  them  own  PC-based  logging  systems. 
The  .ACM  data  presented  here  also  recorded  two 
“slow"  channels  of  data  taken  from  the  UD.ATS 
system  to  ensure  exact  temporal  sychronisation 
between  the  different  logging  systems. 

The  .ACM  was  mounted  60  cm  above  the  sea 
bed  with  its  transducer  axis  vertical  on  the  hori¬ 
zontal  arm  of  an  H-frame  100  m  from  MWL 


(b)  / 


Fig. 2.  (a)  Cross  secaon  of  Ouecnsiand  beach  with  the  four 
major  instrument  locations.  (bi  H-Trame  w-ith  acoustic  concen¬ 
tration  meter  lAC.Mi  and  eiectromacnctic  current  meters 
lEMCSl. 
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Fig. 3.  Sieved  sand  size  distribution  in  half  phi  intervals  of 
surricial  sample  taken  from  Queensland  beach  close  to  the 
■ACM  location. 

(Fig.ZI.  Current  meters  were  fixed  to  one  of  the 
vertical  arms  of  the  frame  at  heights  of  20  cm  and 
50  cm  above  the  bed  so  that  they  measured  the 
current  close  to  the  area  ensonified  by  the  ACM. 

The  instrumentation  was  positioned  by  divers 
who  also  made  visual  observations  of  bottom 
roughness  and  ripple  wavelength  whenever  pos¬ 
sible.  .A  video  camera  was  mounted  near  the  90  m 
frame  and  recordings  made  of  resuspension  events, 
but  no  direct  information  about  ripple  dimensions 
could  be  gathered  for  this  storm  due  to  poor 
visibility  and  low  light  levels. 

3  The  acoustic  measurement  of  suspended 
sediment 

3.1  The  acoitsik  concentration  meter 

The  acoustic  concentration  meter  (ACM)  used 
in  this  field  experiment  is  similar  to  that  described 
by  Huff  and  Fiske  (1980).  Hanes  et  al.  (1988)  and 
Libicki  et  al.  (1989).  The  acoustic  frequency  is  2.8 
,\lHz  and  the  pulse  length  is  about  10  ps.  Range 
gating  the  backscaitered  acoustic  signal  allows  the 
sediment  concentration  profile  to  be  estimated  at 
124  range  "bins",  with  a  vertical  resolution  of 
O.Si  cm.  The  pulse  repetition  rate  was  10  Hz  and 
four  profiles  were  averaged  before  storing  the  data 
on  disk  together  with  two  values  from  "slow" 
channels,  on  this  occasion  two  components  of  an 


electromagnetic  current  meter.  Disk  access  time 
and  other  software  constraints  meant  that  one 
average  profile  was  recorded  every  0.58  s.  1250 
average  profiles  (12  min)  were  recorded  during 
each  run. 

3.2  Calibration  of  the  .ACM 

The  ACM  was  calibrated  in  a  laboratory  recircu¬ 
lation  tank  (Vincent  et  al..  1986)  using  sand  taken 
from  the  beach  at  Queensland.  The  calibration 
curves  for  sand  at  four  concentrations  are  shown 
in  Fig.4a.  the  digital  value  of  the  nominal  voltage 
(8  bit.  0-255  scale)  plotted  against  distance.  The 
concentrations  are  those  measured  by  withdrawing 
a  water  sample  from  the  centre  of  the  recirculation 
tank  and  then  filtering  and  weighing  the  filtrate. 
These  calibration  profiles  are  each  an  average  of 
1250  individual  profiles.  The  acoustic  backscatter 
intensity  V-  (expressed  here  in  terms  of  the 
notional  voltage  V  seen  by  the  system)  from  a 
uniform  field  of  scatterers  of  concentration  C{r)  is 
assumed  to  be  an  inverse  function  of  range  r  with 
corrections  for  attenuation  due  to  water  and  to 
the  scatterers: 

=  Br'-C(ri  e.xp  |-4xr  -  J  ^C(ridr|  (1) 

where  B  is  a  constant  which  includes  the  beam 
strength  and  sediment  backscattering  cross  section, 
a  is  the  attenuation  due  to  water  and  /J  is  the 
attenuation  of  the  beam  due  to  the  presence  of 
sediment.  Here  we  assume  a  value  of  1.1  dB/m  for 
the  absorption  of  2.8  MHz  sound  by  seawater  and 
that  0  is  2.5  X  10"^cm"‘  (mg  I"')"'  (Sheng  and 
Hey.  1988).  B  depends  on  the  acoustic  pow'er 
emitted  by  the  transducer,  the  sediments  in  the 
water  column  and  the  gain  of  the  receiver  circuits. 
For  this  calibration  (using  Queensland  sand)  a 
value  of  6.4  x  (mV-  cm-  mg"  ‘  !)  was  found  to 
fit  the  calibration  data  best.  Normalising  the  back- 
scattered  acoustic  energy  using  the  above  equation 
reduces  the  calibration  data  to  that  shown  in 
Fig.4b:  the  region  close  to  the  transducer  (5-25  cm) 
corresponds  to  the  near-field  region  and  fits  the 
theoretical  extent  of  the  near  field  well  ika~.  where 
k  is  the  acoustic  wavenumber  and  a  is  the  radius 
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Oistance  in  cm.  from  transducar 


Fig, 4,  (a)  Calibration  curves  for  sand  from  Queensland  beach  tn  a  recirculation  lank  for  four  concentrations,  (b)  Normalised 
backscatier  signal  intensities  from  the  four  calibration  curves  above.  The  near-field  extends  to  25  cm.  Dashed  lines  are  ±  10%  of 
the  normalised  value. 


of  the  transducer).  The  field  data  are  converted  to 
concentrations  using  a  step  wise  application  of 
eqn.  (10)  with  a  linear  correction  in  the  near-field 
region  as  indicated  in  Fig.4b. 

3.3  Limitations  of  the  A  CM 

There  are  limitations  to  the  effective  use  of  the 
.ACM  for  measuring  rapid  changes  in  the  concen¬ 
tration  of  sediment  above  the  sea  bed  (Vincent  et 
al.,  1986).  Breaking  waves  inject  bubbles  into  the 
water  column  and.  as  air  bubbles  have  a  high 
acoustic  impedance,  the  .ACM  can  clearly  "see" 
them  but  it  is  difficult  in  any  quantitative  manner 
to  differentiate  between  suspended  sand  and  air 
bubbles;  hence  the  .ACM  should  not  be  used  near 


the  surf  zone.  .At  high  suspended  concentrations, 
difficulties  can  also  be  encountered  because  the 
acoustic  beam  is  attenuated  so  rapidly  that  distance 
with  little  energy  penetrates  to  the  full  range  of 
the  ACM;  at  100  mg  1~*  the  attenuation  is  0.6  db 
m'*  but  provided  high  concentrations  (>  1  g  1"‘) 
are  encountered  only  within  a  few  centimetres  of 
the  sea  bed.  reasonable  estimates  of  the  concentra¬ 
tion  can  still  be  made. 

3.4  Sensitivity  to  particle  size 

Calibration  tank  tests  with  sieved  sands  suggest 
that  an  ACM  operating  at  2.8  MHz  is  most 
sensitive  to  sands  between  1.0  and  1.50  (0.5- 
0.35  mm).  The  calibration  equation  described 
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above  assumes  that  the  size  distribution  of  the 
suspended  sediment  is  always  the  same  as  that  of 
the  surficial  sample  collected.  In  the  field  the  size 
distribution  will  vary  both  spatially  and  tempo¬ 
rally  so  the  calculation  of  the  suspended  concen¬ 
tration  from  the  backscattered  acoustic  intensity- 
may  involve  considerable  uncertainty.  The  magni¬ 
tude  of  this  uncertainty  cannot,  therefore,  be 
precisely  quantified  but  for  sand  sizes  which  are 
much  less  than  the  wavelength  of  the  sound 
(which  is  0.58  mm  at  2.8  MHz);  halving  the  sand 
size  halves  the  intensity  of  the  backscattered 
signal. 

5.5  Organic  material  in  suspension 

A  further  confusion  in  the  interpretation  of  an 
■ACM  record  was  highlighted  during  a  pre-storm 
trial  of  the  instrumentation.  Wave  heights  were 
very  low  (<10cm)  and  divers  had  reported  that 
oniy  occasional  sediment  suspension  was  occurring 
from  the  top  of  ripples  on  the  sea  bed.  Visual 
monitonng  of  the  .ACM  signal  during  this  period 
showed  clear  and  regular  acoustic  returns  from 
the  lowest  5  cm  of  the  water  column,  returns  which 
varied  in  intensity  over  a  wave  period  as  would 
be  expected  from  sediment  resuspension  by  waves. 
A  further  inspection  of  the  sea  bed  beneath  the 
.ACM  by  one  of  the  authors  revealed  that  there 
was  a  layer  of  fine  organic  material  floating  close 
to  the  sea  bed  and  being  advected  through  the 
.ACM  beam  by  the  wave  motion.  Spatial  variations 
in  the  density  and  thickness  of  this  organic  material 
were  seen  by  the  ACM  as  a  temporal  variation  in 
suspended  concentration  equivalent  to  a  sand  sus¬ 
pension  of  about  100  mg  1" '  in  the  lowest  3-5  cm 
of  the  water  column. 

There  is  little  that  can  be  done  to  eliminate  this 
probiem  as  in  many  circumstances  it  will  not  be 
known  if  organic  material  is  present  in  the  water. 
During  the  deployment  described  below  we  expect 
the  organic  material  to  be  dispersed  more  uni¬ 
formly  through  the  water  column  than  in  the 
situation  d-escribed  above.  It  should  also  be  noted 
that  a  uniform  suspension  will  appear  as  a  signal 
in  the  transport  due  to  the  mean  current  but  should 
integrate  to  zero  for  wave-induced  transport, 
assuming  the  particles  move  with  the  water. 


3.6  Accuracy  and  precision  of  the  .ACM 
measurements 

Because  of  the  likely  variation  in  the  sand  size 
distribution  with  height  above  the  sea  bed  it  has 
not  been  possible  to  calculate  the  accuracy  of  the 
suspended  sand  concentrations:  Vincent  and  Green 
(1990)  show  that  the  apparent  concentrations 
above  the  sea  bed  under  waves  decreased  faster 
than  predicted  by  the  Glenn  and  Grant  (1983) 
model  and  suggested  that  this  was  due  to  finer 
material  diffusing  higher  into  the  flow.  During  the 
Queensland  experiment,  suspended  sand  samples 
were  taken  at  a  number  of  heights  above  the 
botlon  near  the  90  m  size  using  a  Nielsen  suction 
sampler:  there  was  some  indication  of  a  decrease 
in  the  model  size  of  suspended  sand  but  the 
differences  were  small  (Greenwood.  1989,  pers. 
commun.).  The  precision  of  the  measurements  is 
believed  to  be  high  and  the  calibration  curves 
suggest  a  precision  of  —  10%  in  the  assessment  of 
the  sand  concentration,  equivalent  to  the  calibra¬ 
tion  sand  size  distribution. 

4  .Acoustic  measurements  of  suspended  sand 

concentrations  during  a  mild  storm  event 

4.1  JA-'ave  and  current  conditions 

The  event  chosen  for  analysis  was  a  25  h  period 
of  onshore  winds  beginning  at  0800  h  on  28th 
October  1987  (Julian  days  301-302).  During  this 
period,  twenty  12  min  runs  of  the  .ACM  were 
recorded,  of  which  sixteen  had  synchronous 
electromagnetic  current  meter  (EMCM)  data;  runs 
are  identified  by  Julian  day  plus  a  sequence  number 
(301  1-301  13.  302 '1-302.7).  Examples  of  the  tem¬ 
poral  and  spatial  variations  in  the  suspended  sand 
concentrations  during  four  runs  (301.3.  301 'll. 
302  1  and  302  5)  are  shown  in  Fig.5  together  with 
the  EMCM  current  speeds.  Mean  currents,  wave- 
induced  currents  and  wave  period  for  the  complete 
period  are  summarised  in  Table  1. 

Speeds  from  the  shore-normal  component  of 
the  E.VICM  at  20  cm  above  the  bed  I  waves  were 
approaching  directly  onto  the  shore)  were  used 
to  calculate  the  wave-induced  current  and  the 
wave  period:  no  significant  difference  was  found 
between  the  speeds  measured  by  the  EMCMS  at 
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20  cm  and  50  cm.  Mean  currents,  speed  and  direc¬ 
tion.  were  calculated  from  the  two  components 
of  the  EMCM  at  20  cm.  The  method  of  analysis 
follows  that  used  by  Vincent  and  Green  (1990). 
After  subtracting  the  mean  current,  the  wave 
period  is  found  from  the  downward  zero-crossing 
period.  The  maximum  current  is  then  found  for 
each  wave  half-cycle  and  a  "significant  current" 
(defined  in  a  similar  way  to  the  significant  wave 


height,  as  the  average  of  the  highest  one-third  of 
the  maximum  currents)  is  identified  for  the 
onshore  half-cycles  and  for  the  offshore  half¬ 
cycles  of  the  w  ave  (Table  1 ). 

•t.2  Position  of  the  sea  bed 

Suspended  sand  concentrations  decrease  rapidly 
away  from  the  bed  so  an  accurate  assessment  of 
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the  position  of  the  sea-bed  is  required.  The  level 
of  ma.ximum  acoustic  backscaiter  was  interpreted 
as  an  echo  from  the  sea  bed;  this  level  remained 
constant  ( -  1  bin.  r  0.8  cm)  throughout  each  run. 
indicating  that  any  bed  features  present  were  not 


moving  significantly  during  that  time,  or  that  the 
amplitude  of  the  bedforms  was  small  (~1  cm). 
The  acoustic  footprint  of  the  ACM  is  about  3.5  cm 
at  the  60  cm  range  (defined  by  the  position  of  the 
-20db  intensity  contour  of  the  acoustic  beam). 


(C) 
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Fic.5  Ore;  -^cale  plots  of  the  suspenoeC  sand  concentrations  img  ll  and  the  current  speed  tin  the  cross-shore  directioni  for  four 
bursts,  la !  .'01  t  (bi  301  11.  icj  302  I .  idi  302  5. 
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TABLE  1 


Wave  and  current  conditions  during  the  24  h  mild  storm  event,  showing  the  mean  currents  and  the 
significant  (sig)  and  average  (av)  wave-induced  maximum  currents  for  the  onshore  and  offshore  half- 
cycies  lincm  s*').  Bursts  for  which  sychronous  current  data  were  not  available  are  shown  with 
astensks 


Burst 

No. 

Time 

(h) 

Mean  currents  (cm  s"') 

Wave-induccQ 

currents 

Wave 

period 

(SI 

On  off 

Along 

Onshore 

Slg 

av 

Offshore 

5ig 

av 

301  1* 

0815 

2.6 

-  1.1 

25 

16 

-27 

-  16 

4.70 

301  2* 

0900 

2.6 

-  1.1 

''n 

16 

_  ■y“t 

-  16 

4.70 

301  3 

1100 

2.3 

-0.6 

24 

15 

-24 

-  15 

4.36 

301  4* 

1300 

3.1 

-1.6 

30 

19 

-31 

-  19 

4.60 

301  5 

1530 

4.1 

-1.9 

39 

26 

-40 

—  25 

4.73 

301  6’ 

1645 

4.1 

-1.9 

39 

26 

-39 

—  25 

4.74 

301  7 

1730 

5.8 

-2.8 

52 

35 

-  S2 

-32 

4.92 

301  8 

1830 

4.8 

-1.6 

58 

39 

-  56 

-35 

4.95 

301  9 

1930 

7.5 

-4.3 

61 

39 

-55 

-34 

4.95 

301  10 

2030 

9.1 

-5.8 

73 

48 

-67 

-40 

5.12 

301  11 

2130 

6.8 

"3.7 

61 

40 

-57 

-35 

5.08 

501  12 

2200 

6.0 

-3.8 

62 

40 

-58 

-37 

5.34 

301  13 

2330 

5.7 

-3.2 

50 

32 

-47 

-30 

5.54 

302  1 

0130 

4.2 

-2.0 

37 

23 

-  37 

-23 

5.10 

302  2 

0230 

2.8 

0.5 

43 

25 

-41 

-23 

5.85 

302,3 

0330 

2.6 

0.7 

37 

24 

-  36 

6.12 

302  4 

0530 

-.3 

1.2 

40 

26 

-39 

-25 

6.30 

302  5 

0630 

2.0 

0.9 

32 

20 

-30 

-  19 

5.91 

302  6 

0730 

1.2 

37 

25 

-33 

-21 

6.40 

o 

J 

0850 

2.9 

1.2 

38 

24 

-34 

-21 

6.38 

A  time  series  of  bed  heights  from  the  individual 
profiles  of  run  301,'8  (Ar=0.58s)  and  profiles 
smoothed  over  ten  pulses  (Ar=5.S  s)  is  shown  in 
Fig.6. 

Over  the  period  of  the  storm  event  however,  the 
run-average  sea  bed  position  varies  over  more  than 
4  cm  (Fig. 7)  but  it  is  not  possible  to  conclude  how 
much  of  this  variation  is  due  to  local  erosion  or 
deposition,  or  to  the  position  of  bedforms  relative 
to  the  acoustic  beam  axis. 

4.3  Run-awra^c  suspended  sand  conceniraiions 

The  concentrations  of  sand  in  suspension  vary 
rapidly  from  second  to  second,  producing  a  picture 
of  complex  suspension  events  (Fig.5).  Time  averag¬ 
ing  the  concentration  profiles  over  the  period  of  a 
run  showed  an  unexpected  pattern  of  suspended 
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Fig.6.  (a)  Time  senes  of  the  posiiior.  of  the  sea  bed.  identified 
by  the  maximum  echo  strength  for  the  first  210  s  of  burst  201  8 
(A7'=0.58sl.  ibi  Time  senes  averaged  over  ten  profiles  iAT= 
5.S  St. 
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Fi".'  Bursi-averagsc  posmons  oC  ;he  seu  btc  over  the  period 
of  the  storm. 

sand  concentrations,  with  the  lowest  concentra¬ 
tions  apparently  associated  with  the  period  of 
highest  waves  (shown  here  by  the  largest  wave- 
induced  current  speeds).  A  time  senes  of  run- 
averaged  concentrations  at  four  heights  (2.5,  5,  12 
and  25  cm)  together  with  the  wave-induced  current 
speed  and  wave  period  is  shown  in  Fig.8.  There  is 
a  negative  correlation,  significant  at  99%.  between 
the  significant  wave  orbital  speed  and  concentra¬ 
tions  at  2.5.  5  and  12  cm. 

Three  average  run-mean  profiles  have  been 
plotted  in  Fig. 9  for  periods  when  profile  shape 
was  not  changing  rapidly:  the  rising  phase  of  the 
storm  (301  1-301  4).  the  maximum  wave  height 
phase  (301  "-30112)  and  the  waning  phase 
(302  2-302  ").  These  average  profiles  differ  from 


's  '2  '6  20  24  4  e 


0?  30'  23  3C2 

Fil  S  f  j  I  Burst- jicraged  suspended  sand  concentrations  at  2.5. 
5.  12  and  25  cm  above  the  sea  bed  (note  the  log  scale),  (b) 
Sicnincant  ivasc-induced  currents  (solid  line)  and  the  wave 
periods  idashed  line;,  ici  ripple  steepness  from  Carstens  et  al¬ 
ii  %4) 


Fig.9.  Average  profiles  for  301.1-301/3  (the  rising  storm). 
301  6-301  II  (maximum  wave  heightsi  and  302.2-302/7  (the 
waning  storm  |. 

each  other  at  most  heights  but  the  important 
differences  are  in  the  lowest  20  cm;  above  this 
level,  much  of  the  backscatter  may  be  due  to 
suspended  organic  material  (see  later). 

5  Wave-current  interaction,  bed  roughness  and  the 

resuspension  of  sand 

Under  a  combination  of  a  steady  current  and 
oscillatory  flows,  which  are  generated  near  the  sea 
bed  by  surface  waves,  an  interaction  occurs 
between  the  waves  and  the  currents.  Close  to  the 
sea  bed  there  is  a  shallow  wave  boundary  layer 
while  above  this  the  steady  current  experiences  an 
apparent  roughness  which  is  enhanced  by  the 
presence  of  the  waves  (Grant  and  Madsen,  1979). 
We  have  used  the  model  developed  by  Grant  and 
Madsen  (1979.  1982)  and  Glenn  and  Grant  (1983) 
for  the  wave-current  interaction  over  a  mobile 
bed,  including  the  suspension  of  sand  and  a  correc¬ 
tion  for  anv  sediment-induced  self-stratification 
which  mav  occur  (this  mode!  is  henceforth  referred 
to  as  the  Grant-.Madsen-Glenn  or  GMG  model), 
to  predict  the  combined  flow  characteristics  and 
the  bedforms.  We  have  also  compared  the  concen¬ 
tration  measured  2  cm  above  the  sea  bed  with  that 
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predicted  by  the  model.  The  model  integrates  the 
Smith  and  McLean  (1977)  expression  for  the  refer¬ 
ence  concentration  Cir,,)  over  a  wave  period  T: 


,”r  =  T 

.1  =  0 


1  -  ■,  nS'i:i 


dr 


(2) 


where  C,,  is  the  bed  volume  concentration  of 
sediment.  Vq  is  an  empirical  constant.  Cq  is  the 
roughness  height  in  the  wa\e  boundary  layer  and 
S'ln  is  the  instantaneous  normalised  excess  skin 
friction,  de.hned  as: 


S'(f) 


Ur  is  the  sand  fall  velocity  and  •;  is  an  empirical 
constant  set  to  0.74  following  Businger  et  al's. 
(1971)  work  in  the  atmospheric  boundary  layer. 

The  ripple  geometry  is  defined  empirically  from 
the  laboratory  data  of  Carstens  et  al.  (1969).  The 
ratio  of  the  maximum  Shields  parameter  tj/‘  over 
a  wave  cycle  to  the  critical  Shields  parameter  tj/^ 
is  used  to  define  the  break-off  point  where  the 
npple  steepness  begins  to  decrease: 

(w' l.S  S.- "  (10) 

5,  is  a  dimensionless  sediment  parameter  given 
by: 

S.  =  (Z»4v)[(5-l)gZ)f-*  (11) 


where  :'(/)  is  the  instantaneous  skin  friction  and 
IS  the  critical  shear  stress  for  initiation  of 
motion.  Various  authors  have  suggested  different 
values  for  (Dyer.  1980:  W’iberg  and  Smith.  1983; 
Hill  et  al..  1988:  Vincent  and  Green.  1990)  and 
recently  Drake  and  Cacchione  (1989)  presented 
data  which  showed  that  -/q  decreased  rapidly  as 
excess  skin  friction  increased. 

The  steady  current  profile  for  non-straiified  con¬ 
ditions  appiicabie  to  all  the  data  presented  here, 
both  within  (r<()*)  and  above  the  wave 

boundarx-  layer,  is  given  by: 


:  ln(--  Co).  -<5. 

(4) 

K')  Inir.CoJ.  Od, 

(5) 

where  r,,  =  A(,  .’0  and  is  the  hydrodynamic 
roughness,  and  Coc  is  the  apparent  roughness  that 
retards  the  mean  flow  above  the  wave  boundary 
layer.  The  metion  velociix  characteristic  of  c>c)„ 
is  (.  .c  =  (ir:.  p)"  '.  where  is  the  time-averaged 
instantaneous  shear  stress  at  the  top  of  the  bound¬ 
ary  layer.  =  (‘T-»i  ff'  is  the  friction  velocity 
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in 

lane 

ous  boundary  shear  stress. 

T 
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trai  concentration  profiie  is: 
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1  = 

C;-'. 

1C  c,,)''.  c<J„ 
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1  = 

Cl 
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y.  ih 

the  suspension  parameter: 

y.  =  ■ 

u . 

for  c<d„ 

(8) 

a  =  ■ 

■u , 

,  kC 

o' 

•1 

S' 

(9) 

where  v  is  the  kinematic  viscosity,  s  is  the  specific 
gravity  of  sand,  g  is  the  acceleration  due  to  gravity 
and  D  is  the  sediment  diameter.  The  Shields  param¬ 
eter  is  defined  in  the  usual  way  as: 

ij/  =  r:pis- \}"D  (12) 

where  r  is  the  appropriate  bed  shear  stress  (the 
skin  friction)  and  p  is  the  density  of  seawater.  The 
ripple  heights » (from  table  1  in  Grant  and  Madsen. 
1982)  are  gi\en  by: 

t]  a  ==0.22  (v  for  (13) 

);;a  =  0.48  \u/'  for 

(14) 

where  a  is  the  wave  orbital  excursion  amplitude. 
The  ripple  lengths  /.  are: 

;..'^  =  6.25  (y'  1^,)°  °“  for  {ijj'  wJb  (15) 

/.  f?  =  3.6  S*”'’-'’  (li/'  for  (l^''l/'J>(l^'  ’l/'c)B  (16) 

We  have  run  the  GMG  model  using  the  wave  and 
current  data  from  the  EMCMs;  the  onshore  signi¬ 
ficant  wave-induced  current  was  used  to  represent 
the  wave  conditions  in  each  run.  The  model  was 
allowed  to  denne  its  own  equilibnum  bedform 
conditions  and  to  correct  for  any  self-stratification 
(no  significant  sediment  stratification  was  encoun¬ 
tered  during  any  of  the  runs  presented  here  and 
ail  the  results  are  effectively  for  neutral  conditions). 
Table  2  shows  the  npple  height  rj.  wavelength  A 
and  steepness  A.  the  values  of  L',cw  -oc  ^^d 
Co  and  the  excess  skin  faction  predicted  by  the 
model,  the  last  term  based  on  the  assumption  that 
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TABLE  2 


Predicted  bedforms.  fnction  velocities  and  roughness  lengths  from  the  Grant-Madsen-Glenn  model  for 
each  run.  and  the  normalised  excess  bed  friction 


Burst 

No. 

Bedforms 

(cm) 

icm  s) 

(cm  s) 

-'»C 

icm) 

-O 

(cm) 

Normalised 
excess  skin 
friction 

'0 

n 

(cm) 

(cm) 

n  '■ 

301  1 

3.3 

21 

0.16 

4.2 

0.64 

8.0 

3.6 

0.48 

301  2 

3.3 

21 

0.16 

4.2 

0.64 

8.0 

3.6 

0.48 

301.3 

3.4 

0.15 

4.4 

0.57 

8.0 

3.8 

0.51 

0.5 

9.7 

301  4 

3.9 

0.14 

5,4 

0.91 

9.3 

4.3 

0.54 

301  '5 

2.6 

28 

0.093 

5.2 

1.05 

8.6 

3.6 

0.25 

2.4 

4.6 

301,6 

2.6 

28 

0.093 

5.2 

1.05 

8.6 

3.6 

0.25 

3017 

1.7 

31 

0.056 

5.6 

1.38 

9.0 

3.2 

0.14 

4.6 

0.72 

301:8 

1.5 

31 

0.047 

6.0 

1.21 

9.5 

3.7 

0.13 

5.7 

0.50 

301,9 

1.3 

31 

0.042 

6.2 

1.82 

9.9 

3.1 

0.13 

6.5 

0.35 

301'I0 

1.0 

33 

0.030 

7.6 

2.42 

11.7 

3.4 

0.14 

7.6 

0.20 

301,11 

1.4 

32 

0.042 

6.4 

1.69 

9.9 

3.3 

0.13 

6.4 

0.17 

301.12 

1.4 

34 

0.042 

6.5 

1.64 

9.9 

3.5 

0.13 

6.5 

0.27 

301,13 

2.1 

35 

0.06  i 

6.3 

1.47 

8.9 

3.4 

0.17 

4.1 

0.46 

302! 

3.2 

31 

0  'C.' 

5.7 

1.13 

8.8 

3.9 

0.32 

2.0 

0.29 

302.'2 

3.1 

37 

lOO 

6.6 

0.80 

8.8 

4.9 

0.26 

2.8 

1.70 

302,'3 

4.2 

38 

O.IIO 

7.1 

0.84 

9.1 

5.5 

0.44 

1.9 

1.63 

302,4 

3.8 

0.097 

7  ^ 

0.81 

9.0 

5.5 

0.36 

2.3 

1.95 

302,'5 

5.0 

36 

0.140 

7  '* 

0.73 

9.6 

6.0 

0.67 

1.2 

1.54 

302,'6 

4,4 

40 

O.IIO 

7.5 

0.91 

9.2 

5.* 

0.47 

1.8 

2.16 

302,7 

4.2 

40 

0.105 

7.4 

0.95 

9.1 

<  < 

0.43 

2.0 

1.84 

the  average  skin  friction  over  ripples  is  the  same 
as  for  that  over  a  flat  bed  (Grant  and  Madsen, 
1982). 

5.1  Variations  in 

The  values  of  Vo  shown  in  Table  2  are  the  values 
required  by  the  model  to  match  the  predicted 
suspended  sand  concentration  at  2  cm  above  the 
sea  bed  with  the  measured  values  (average  of 
concentrations  from  bins  2  and  ?  above  the  bed). 

reference  concentration  at  2  cm.  is  preferred  to 
that  at  !  cm  because  of  the  uncertainty  of  r  1  bin 
(  rb-S  cm  I  in  the  bed  position;  however,  the  strong 
concentration  gradient  near  the  bed  will  mean  that 
C(Zo)  is  very  sensitive  to  the  matching  height 
chosen.  The  variation  in  Vq  with  e.\cess  skin  friction 
is  shown  in  Fig.  10;  the  trend  is  similar  to  that 
found  by  Drake  and  Cacchione  (1989)  but  the 
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Fig.  10.  Plot  of  '.he  variation  in  the  resuspension  coelhcienl  Vp 
against  the  ma.ximum  (over  one  wave  period)  excess  skin 
friction.  •This  study;  C Vincent  and  Green  il990);  -Drake 
and  Caccione  ( 1989). 
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values  of  Vo  are  generally  an  order  of  magnitude 
higher.  Drake  and  Cacchione's  values,  however, 
were  extrapolated  from  nephelometer  data  taken 
at  a  height  of  1.9  m  over  a  silt  bed  with  a  mean 
grain  size  of  16-25  pm.  Towards  the  end  of  the 
storm  period  when  wave  heights  were  decreasing 
and  bed  roughness  was  greatest,  the  values  of  Vj, 
were  close  to  the  value  of  2  x  10”^  found  by  Smith 
and  McLean  (1977)  in  the  Columbia  River.  During 
the  higher  excess  skin  friction  period  the  values 
were  ver\'  close  to  those  measured  by  Vincent  and 
Green  (1990)  for  resuspension  on  a  beach  with 
similar-sized  sand. 

5.2  Comparison  of  modelled  and  observed 
concentration  profiles 

There  are  important  differences,  similar  to  those 
noted  by  Vincent  and  Green  (1990).  between  the 
average  concentration  profiles  and  the  concentra¬ 
tion  profiles  predicted  by  the  GMG  model  (Fig.  1 1 ). 
In  particular,  the  observed  gradients  of  suspended 
sand  are  steepest  close  to  the  bed  (in  the  wave 
boundary  layer),  whereas  the  model  predicts  a 
much  more  uniform  gradient  close  to  the  bed. 
.Also,  it  is  difficuilt  to  identify  any  position  as  the 
top  of  the  wave  boundar}'  layer  in  the  observed 
profiles. 

While  these  fundamental  differences  in  the  pro¬ 
file  shapes  show  that  the  GMG  model  does  not 
satisfactorily  predict  the  average  naturally  occur¬ 
ring  sand  concentration  profiles  in  this  nearshore 
environment,  the  model’s  shortcom.ings  are  not 


unexpected.  The  GMG  model  is  a  two-layer  model 
with  the  vertical  eddy  diffusivity  being  time-invari¬ 
ant  and  linear  with  height:  it  uses  a  single  wave 
height,  wave  period,  ripple  height  and  ripple 
length.  .A  real  sea  has  a  spectrum  of  wave  heights 
and  periods  so  the  height  of  the  wave  boundary 
layer  will  var\-  from  wave  to  wave,  the  ripples  too 
are  likely  to  be  responding  dynamically  to  the 
changing  wave  characteristics  and  vortex  ejection 
per  se  cannot  be  considered  because  of  the  form 
of  the  eddy  diffusivity  profile.  The  differences  do 
nevertheless  warn  of  problems  in  using  the  GMG 
model  to  predict  suspended  sand  concentration 
profiles. 

6  On-offshore  suspended  transport  rates 

The  transport  of  sand  can  be  divided  into  two 
components,  an  adveclive  term  and  a  wave  oscilla¬ 
tory  term.  The  advective  component  of  transport 
is  due  to  the  steady  (time-averaged)  current  U(:) 
and  the  entire  suspension  of  sand  in  the  water  is 
bodily  transponed.  In  the  wave  oscillatory  compo¬ 
nent,  where  the  net  transport  of  water  is  by 
definition  zero,  sand  transport  occurs  due  to  time 
lags  between  fiuctuations  in  the  current  U„{z,t) 
and  fluctuations  in  suspended  sand  concentrations 
C(c,/).  While  the  instantaneous  current  can  be 
divided  easily  into  steady  and  wave  oscillatory 
components.  Hz.t)=U^(z,t)+U{z),  it  should  be 
noted  that  the  same  cannot  be  done  for  the  sus¬ 
pended  sand  concentration  in  any  meaningful  way: 
the  sand  concentration  at  any  instant  is  due  to  the 


Fic. II.  Compansor  betwcsn  the  measurec  barst-averaged  concentration  profiles  (solid  lines)  and  those  predicted  using  the  GMG 
model  with  nppie  dimensions  determinea  from  Carsiens  et  al.  (1969)  (dashed  lines)  and  from  Nielsen  (1981)  (dotted  lines)  for  bursts 


CONTROL  OF  SUSPENDED  SAND  CONCENTRATION  BY  BED  ROUGHNESS 


13 


excess  skin  friction  due  to  both  waves  and  currents 
(alongshore  as  well  as  on-offshorej  over  a  mobile 
bed  formed  by  the  combined  system.  Nevertheless, 
it  is  useful  to  compute  the  relative  importance  of 
the  two  components,  as  defined  below. 

The  on-onshore  component  of  the  transport  of 
suspended  sand,  averaged  over  the  period  of  a 
burst,  has  been  calculated  in  several  ways,  as 
follows: 

(1)  For  steady  currents,  as  a  function  of  height 

above  the  bed  and  integrated  over  the 

water  column 

0(.-|„„„  =  r(.-i  G--)  (17) 

and 

Q„cun=  "l""  Q(-Ua„  (18) 

z  =  >ca  oed 

w  here  Gci  is  the  burst-mean  concentration  profile. 

(2)  For  the  wave  oscillatory  currents,  again  as 
a  function  of  height  from  the  bed  (2(-)wave)  and 
integrated  over  the  water  column  (Qwave): 

r*  Hurst  _ 

1  .-iCl-  r)  (19) 

:  -  u 

and 

2wav.=  ’  l'"  G-Uv,  (20) 

r  =  <e3  n«d 

The  average  %elocity  profile  was  that  predicted 
by  the  GMG  model  with  neutral  stratification 
matched  with  the  measured  currents  at  20  cm.  This 
also  assumes  that  the  current  measured  at  20  cm 
is  in  the  potential  flow  region.  Closer  to  the  bed 
in  the  wa\e  boundary  layer,  the  current  velocity 
at  height  c  is  given  by: 

ker  2,  c  i  kei  2,  c  \ 

:i  =  I  1  -  .  .  : 

Co  *  I  kei  Co,/ 

>.  C'..ix.ri  (21) 

where  x./i  is  the  wave-induced  current  outside 
the  w  ave  boundary  layer.  c„  =  Cq  ( tind  c  =  z.l  where 
n.  and  kcr  and  A'c/are  Keh'in  functions 
of  the  first  order  (Grant  and  Madsen.  1978).  This 
equation  snows  the  phase  of  the  current  changing 
towards  the  sea  bed:  the  magnitude  of  this  change 
Is  about  20  at  1  cm.  less  than  the  time  resolution 


of  the  measurements  (0.58  Hz)  for  most  waves,  so 
no  attempt  was  made  to  correct  for  phase  shift  in 
these  results. 

6.1  Veriically  integrated  transport 

The  vertically  integrated  values  of  Smean  ^nd 
2w»ve^  which  describe  the  total  sand  transport 
(flux)  past  the  measurement  position,  divided  into 
the  two  components  described  above,  are  shown 
in  Fig.  12.  The  steady  currents  have  an  offshore 
component  throughout  the  experimental  period  so 
the  transport  due  to  the  steady  current  2me»n  'S 
offshore  throughout.  The  transport  due  to  the 
wave  oscillatory  term  2wave  'S  of  a  similar  magni¬ 
tude  to  2mcan  ^ut  is  generally  directed  towards  the 
shore;  the  uncertainty  limits  shown  are  those  which 
would  result  from  an  error  of  r  1  bin  in  selecting 
the  position  of  the  sea  bed.  There  is  considerable 
variability  in  the  transport  flux  from  run  to  run 
although  the  general  trend  is  for  the  fiux  to  become 
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Fig.  15.  Verticaily  imegraiec  values  of  the  '.-oss-shore  sand 
transport  due  to  the  steady  current  and  to  the  waves 

using  the  GMG  model  with  npple  dimensions  from 
Carsiens  et  a!.  (196^1  i#!  and  from  Nieisen  (1981)  (G).  The 
range  of  Q„„  shown  on  the  Carstens  values  are  those  for  an 
uncertainty  of  -  1  b:r.. 
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more  negatne  (shorewardsj  with  :he  larger  waves. 
The  variability  results  partly  from  the  flux  being 
the  resultant  of  large  transports  in  opposite  direc¬ 
tions  (gross  transport  is  an  order  of  magnitude 
larger  than  the  net  transport!  but  also  from  the 
groupmess  of  the  wavest  the  influence  of  groupi- 
ness  can  be  seen  in  Table  3  where  the  on-offshore 
flux  has  been  calculated  for  each  third  of  run 
30:  2. 

6.2  Tnvispon  profiles 

For  each  run  the  vertical  profiles  of  suspended 
sand  and  the  total  transport 

Q(r),o,^j  in  the  on-offshore  direction  have  been 
computed.  Marked  changes  occur  during  the 
course  of  the  storm,  particularly  in  the  heights  at 
which  maximum  wave  transport  Olcj^ave  ob¬ 
served.  With  the  exception  of  the  last  three  runs 
(302  -i-302  7).  very  little  wave-induced  transport 
occurs  above  20  cm:  the  steady  transport  above 
20  cm  observed  in  a  number  of  runs  is  probably 
due  to  contamination  of  the  backscaitered  signal 
b>  suspended  organic  material.  Figure  13  shows  a 
sequence  of  suspended  transport  profiles  from 
ever}'  other  run  for  the  lowest  20  cm.  The  wave 
component  Qi-fwavt  (cross-hatched  for  clarity)  is 
shorewards  for  all  except  the  first  profile  (run 
301  3).  and  has  a  peak  transport  elevation  which 
decreases  with  height  as  the  waxes  increase  in 
height.  During  the  final  four  runs  i302  4-302.  7) 
transport  occurs  to  a  much  greater  height  from 
the  bed  i~50cm).  although  still  decreasing  with 


height:  this  max  be  due  to  the  efficient  ejection  of 
vortices  from  the  ripple  crests  (predicted  to  be 
steep  at  this  time)  or  perhaps  to  patchy  organic 
material  in  suspension  (see  Fig.cd). 

6.3  Sensitivity  of  the  transport  estimates 

The  sand  transport  values  computed  here  clearly 
depend  on  the  form  of  the  current  profile  (both 
wave-induced  and  the  steady  component)  chosen 
to  occur  beioxx  the  EMCM  at  20  cm.  .All  the  above 
results  use  the  GMG  model  which  includes  ripple 
dimensions  based  mainly  on  laboratory  experi¬ 
ments  by  Carstens  et  al.  (1969).  Howex'er.  Nielsen 
(1981)  noted  that  natural  ripples,  while  obeying 
essentially  the  same  rules  as  laboratory  ripples, 
were  shorter  and  less  steep  due  to  the  irregularity 
of  natural  xvaxes  (Fig.  14).  To  investigate  the  sensi¬ 
tivity  of  the  results  presented  hitherto  to  the  exact 
form  of  the  current  below  20  cm  we  have  tried  a 
number  of  forms  for  the  profile,  based  on  both 
the  GMG  mode!  generating  its  oxvn  ripple  dimen¬ 
sions  (Carstens  et  al..  1969)  and  also  on  the  mode! 
run  with  nppies  defined  from  Nielsen  (1981)  for 
irregular  waves  alone; 

693  -  0.37  in».\/ 

''-a-e.Xp  jQQQ_Q  -5 

»/  a  =  21 

xvhere  a  is  the  xvave  semi-excursion  and  .W  is  a 
mobility  number  defined  as 

.VI  =  («w)‘  (5—  I  ig£> 


TABLE  .t 
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more  negaine  ishorewards)  wuh  ihe  larger  waves. 
The  variabiiitv  results  partly  from  the  flux  being 
the  resultani  of  large  transports  in  opposite  direc¬ 
tions  I  gross  transport  is  an  order  of  magnitude 
larger  than  :he  net  transport)  but  also  from  the 
gri)upiness  of  the  waves;  the  influence  of  groupi- 
ness  can  be  seen  in  Table  3  where  the  on-offshore 
flu.x  has  been  calculated  for  each  third  of  run 
30:  2. 

6.2  Tfanspun  profiles 

For  each  run  the  vertical  profiles  of  suspended 
sand  and  the  total  transport 

'he  on-offshore  direction  have  been 
computed.  Marked  changes  occur  dunng  the 
course  of  the  storm,  particularly  in  the  heights  at 
which  maximcam  vvave  transport  0(rj«.avc  is  ob¬ 
served.  VVith  the  exception  of  the  last  three  runs 
1 302  -i-302  very  little  wave-induced  transport 
occurs  above  20  cm:  the  steady  transport  above 
20  cm  observed  in  a  number  of  runs  is  probably 
due  to  contam.ination  of  the  backscattered  signal 
bv  suspended  organic  material.  Figure  13  shows  a 
sequence  of  suspended  transport  profiles  from 
everv  other  run  for  the  lowest  20  cm.  The  wave 
component  (cross-hatched  for  clanty)  is 

shorewards  for  all  except  the  first  profile  (run 
301  3).  and  has  a  peak  transport  elevation  which 
decreases  with  height  as  the  waves  increase  in 
height.  During  the  final  four  runs  (302  4-302  ”» 
transport  occurs  to  a  much  greater  height  from 
the  bed  ('-50  cm).  although  still  decreasing  with 


height:  this  mav  be  due  to  the  efficient  ejection  of 
vortices  from  the  ripple  crests  (predicted  to  be 
steep  at  this  time)  or  perhaps  to  patchy  organic 
material  in  suspension  (see  Fig.5d). 

6. 3  Sensttiviiy  of  i he  transport  estimates 

The  sand  transport  values  computed  here  clearly 
depend  on  the  form  of  the  current  profile  (both 
wave-induced  and  the  stead.v  component)  chosen 
to  occur  below  -he  EMCM  at  20  cm.  .All  the  above 
results  use  the  GMG  model  which  includes  rippie 
dimensions  based  mainly  on  laboratory  experi¬ 
ments  bv  Carstens  et  al.  (1969).  However.  Nielsen 
(1981)  noted  that  natural  ripples,  while  obeying 
essentially  the  same  rules  as  laboratory  rippies. 
were  shorter  and  less  steep  due  to  the  irregularity 
of  natural  waves  (Fig.  14).  To  investigate  the  sensi¬ 
tivity  of  the  results  presented  hitherto  to  the  exact 
form  of  the  current  below  20  cm  we  have  tried  a 
number  of  forms  for  the  profile,  based  on  both 
the  GMG  model  generating  its  own  ripple  dimen¬ 
sions  (Carstens  et  al..  1969)  and  also  on  the  model 
run  vvith  rippies  defined  from  Nielsen  (1981)  for 
irregular  waves  alone: 

693  -  0.37  ln«  .V/ 
exp  ^  000  -  0.75  In  ’  .Vf 

1/0  =  21 

where  a  is  the  wave  semi-excursion  and  M  is  a 
mobiiitv  number  defined  as 

.V/  =  (owr  (s-  1  igD 
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-  Wave  IransDort 

- Current  transaorl 

-  Total  transDort 


Transoori  rate  tmo'crr  3 


Flu.l?.  Run-avcragcd  suspended  .sand  transoort  profiles  for  the 
lowest  do  cm  ievcr>  other  rum  showir.g  '.he  components  due 
to  t.’ie  steady  current  ana  to  wave.s.  The  cross-hatched  region 
iS  the  shorewarc  transport  due  to  wase;  while  the  diagonal 
hatching  is  the  otTshorc  transport. 

uherd  a;  IS  ihe  wave  radiaa  I'recaency.  The  irans- 
pon  caiculations  have  been  repealed  for  ihe  four 
ume  senes  shown  m  Fie..'  usine  a  '.ariety  of  Sow 
proSies  and  the  results  are  shown  in  Table  4  for 
the  erosshore  component  of  the  transport.  Firstly, 
the  GMG  mode!  was  run  with: 

i  I !  self-adiusting  hed  lopogranhy  lCar^:e^.s  ei  a!..  IVfid). 


302/5 


2Str  Oct  1303011  29th  1JD302) 


Fig.  14.  Companson  of  the  ripple  height  and  length  used  by 
the  GMG  modet  fcr  combined  flow  (based  on  Carsiens  et  al.. 
1969).  and  those  predicted  by  Nielsen  (1981)  for  naturally 
occurring  npp.es  under  waves  alone. 

TABLE  4 

The  sensitivity  of  the  ••ertically  integrated  burst-averaged  trans¬ 
port  computahor.i  to  the  assumptions  about  the  form  of  the 
current  proiiie.  Tne  variants  (I)  and  (2)  and  tav  to  le)  are 
denned  in  the  te.";;.  Transport  values  are  in  mg  cm  s.  For  (c) 
to  teu  shorewarc  transport  is  positive 


lai  (b)  (c)  id)  le) 


301  3 

(1) 

S.5 

20.5 

-S.2 

-5.S 

-  II.O 

(2) 

9.5 

21.6 

-9.7 

-7.1 

-11.0 

301  11 

(1) 

13.0 

20.5 

13.2 

5.2 

15.5 

(2) 

10.9 

17.9 

13.7 

5.8 

15.5 

302  T 

(1) 

!6.6 

43.9 

24.2 

17.6 

28.0 

(2) 

22.1 

52.1 

26.7 

17.6 

28.0 

302  5 

(1) 

3.5 

9.3 

-5.3 

-3.0 

-  10.8 

(21 

-1." 

11. 1 

-9.8 

-4.7 

-  10.8 

(2)  GMG  mode;  with  rippie  dimensions  from  Nieisen  (1981) 
using  wave  orbital  data  only  (no  currents!  and  then  transport 
calculations  made  for  the  steady  current  using 

lai  the  two-i3>er  GMG  flow  profile,  and 
tbi  a  one-iayer  logarithmic  profile  defined  by  the  measured 
current  at  20  cm  and  by  the  npple  height  i;. 

The  wave-mcuced  transport  depends  on  the  profile  in 
the  wave  bouncar;.  layer.  The  following  vananis  were  used: 

(c)  current  profile  denned  as  m  eqn.d  1). 
id)  wave-induced  current  decreases  linearly  from  the  top  of 
the  wave  boundan  layer  to  zero  at  the  bed.  and 
(el  wave-incuced  current  constant  and  equal  to  the  velocity 
in  the  potential  flow  region. 
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Despite  the  large  differences  in  the  nppie  dimen¬ 
sions  predicted  by  the  Carstens  et  al.  ( 1969)  labora¬ 
tory  data  and  Nielsen's  (1981)  held  data,  the 
difference  between  the  transport  rates  (Tabled  (1) 
and  1 2)1  is  tvpicalh  only  20“ o.  The  transport  due 
the  mean  current  using  the  GMG  two-layer 
mode;  gives  values  which  are  about  half  those 
w  hich  are  obtained  if  wave-current  interaction  is 
Ignored  and  a  simple  logarithmic  profile  assumed 
iTabie  lai  and  ibi).  With  the  highest  sand  concen¬ 
trations  occurring  close  to  the  bed  the  form  of  the 
wa\e-!nduced  current  profile  close  to  the  bed 
affects  (2wa%c  significantly  and  is  sensitive  to  the 
form  wave  profile.  The  profile  iTabied  (c))  based 
on  the  Grant  and  Madsen  (1979)  profile  (eqn.ll) 
is  about  mid-way  between  the  resuits  obtained  by 
the  linear  model  id)  and  the  constant  model  (e). 

7  Discussion 

No  direct  measurements  or  observations  of  the 
dimensions  of  the  bedforms  were  available  for  this 
storm  so  it  was  necessary  to  estimate  the  amplitude 
and  wavelength  indirectly.  While  this  is  a  severe 
limitation  of  the  interpretation  of  the  suspended 
sand  data  presented  here  there  is  sumcient  consen¬ 
sus  in  the  literature  to  suggest  that,  at  wave 
conditions  encountered  during  this  storm,  bedform 
steepness  will  decrease  as  wave  heights  increase. 

Following  the  simple  model  for  vortex  shedding 
(Sleath.  1982)  the  vortex  shed  from  the  crest  of  a 
ripple  by  one  wave  half-cycle  is  swept  upwards  by 
the  next  (re\ersed)  half-cycle.  .As  the  ripples 
decrease  in  steepness,  this  ejection  becomes  less 
efficient  and  sand  is  not  lifted  so  high  above  the 
bed.  .At  the  final  stage  of  an  upper  phase  fiat  bed. 
no  vortices  will  'ne  formed  and  all  the  turbulence 
will  be  generated  by  the  oscillatory  fiow  over  the 
sand  grains  resulting  in  a  suspension  remaining 
Closer  to  the  bed. 

Puramciensaium  of  the  rcfcrc/Kc 
conccnirciuon 

Quantitatively  the  changes  in  the  reference  con¬ 
centration  2  cm,  above  the  bed  are  parameterised 
through  the  "constant"  v^.  but  this  is  clearly 
inadccaatc  for  explaining  the  concentration  obser¬ 


vations  shown  in  Fig.S.  The  average  skin  friction 
over  a  rippled  bed  in  the  Grani-Madsen-Glenn 
model  is  assumed  to  be  equal  to  that  which  would 
occur  over  a  fiat  bed  with  the  same  grain  size:  the 
variation  in  the  total  drag  is  assumed  to  be  mainly 
due  to  changes  in  the  form  drag.  These  data 
indicate  that  there  must  be  considerable  variation 
in  the  skin  friction,  averaged  over  a  ripple,  com¬ 
pared  with  a  fiat  bed.  They  are  consistent  with  a 
high  skin  friction  over  the  crests  of  ripples  when 
the  ripples  are  steep,  before  the  break-off  point 
(Grant  and  Madsen,  1982)  is  reached.  Beyond 
break-off  point,  npples  decrease  in  steepness  and 
the  influence  of  ripple  geometry  in  enhancing  skin 
friction  over  the  ripple  crest  decreases  (such  that 
the  ripple-average  skin  friction  also  decreases)  and 
less  resuspension  occurs. 

.A  best-fit  line  t'nrough  the  data  shown  in  Fig.  10 
gives  the  following  relationship  between  Vo  and  the 
maximum  excess  shear  stress  in  n  wave  period; 

Vo  =  7.5x 

It  is  also  noted  that  the  values  of  -/p  in  Fig.  10 
show  a  suggestion  of  hysteresis,  the  higher  values 
of  Vo  occurring  when  the  waves  are  increasing  in 
size  and  hence  the  ripples  would  tend  to  be  steeper 
than  their  equilibrium  value  for  the  waves  and 
sand  suspension  will  be  more  than  expected.  As 
the  waves  decrease  in  height  the  reverse  occurs, 
w'ith  ripples  being  iower  and  hence  less  sand  would 
be  suspended. 

7.2  Transpori  pronies 

Symmetneal  waves  over  symmetrical  ripples  will 
result  in  no  net  sand  transport.  The  transport 
profiles  computed  here,  with  their  maximum  wave- 
induced  transport  between  4  and  10  cm  above  the 
bed.  result  from  fiow  (and  hence  bedform)  asym¬ 
metry.  The  asymmetry  of  the  wave  motion  alone 
is  small,  the  peak  currents  below  the  wave  crest 
being  only  a  few  centimetres  per  second  faster  than 
below  the  trough  (Table  Ij;  however,  from  con¬ 
sideration  of  waxes  alone  the  peak  shear  stress 
(and  hence  peak  sediment  resuspension)  would 
occur  beneath  the  wave  crest.  At  Queensland  beach 
the  mean  on -offshore  current  is  of  the  order  of 
5cm  s’'  and  is  directed  offshore  (Table  1);  this 
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current  and  the  shear  stress  associated  with  it 
appears  be  sumcient  to  bias  the  peak  bed  stress 
towards  the  offshore  wave  half-cycie  (although  the 
shape  of  the  mean  current  prohle  close  to  the  bed 
is  uncertain).  With  such  a  stress  distribution,  the 
time  of  maximum  resuspensior.  would  occur  under 
the  wave  trough  (correlating  with  the  tendency 
seen  in  Fig.}  for  transport  closest  to  tne  bed  to  be 
offshore  I  and  ripples  would  have  tneir  a.symmetry 
directed  offshore.  The  shoreward  .dux  from  5  to 
15  cm  IS  cue  to  the  half-cycle  following  the  maxi¬ 
mum.  .^tre.-.'  sweeping  the  'uspended  sediment  cloud 
uowards  into  the  flow.  The  sand  fall  velocity  of 
2..'  cTi ‘  and  wave  period  o*'  5-6  s  are  consistent 
with  m.aximum  transport  at  the  f-lOcm  height. 
The  importance  of  the  mean  current  to  the  direc¬ 
tion  of  the  wave-induced  transpon  is  also  consis¬ 
tent  with  the  transport  profiles  obserx'ed  by  Vincent 
and  Green  (1990)  wmere  the  mean  flow  was  shore¬ 
ward  and  the  wave-induced  transport  pattern  was 
reversed  (onshore  close  to  the  bed  and  offshore 
between  5  and  15  cm). 

“.c  Time  and  space  structure  of  suspension  and 
~>ed  rougiiness 

The  influence  of  the  groupiness  of  waves  on 
sediment  resuspension  has  been  noted  by  several 
authors  (Clarke  et  al..  1982;  Vincent  et  al..  1982) 
and  can  be  seen  in  some  of  the  suspension  records 
(Fig.f  i.  The  variation  in  ripple  geometry  may  help 
to  explain  this  phenomenon.  Dunng  the  passage 
of  a  group  of  waves  the  sea  bed  will  first  be  under 
the  influence  of  small  waves  (generating  steep 
ripples  at  equilibrium)  and  then  under  larger  waves 
when  ripples  would  be  less  steep  (assuming  the 
break-off  point  had  been  exceeded).  .^  time  lag 
would  be  expected  in  the  response  of  the  bed  to 
the  waves,  so  the  big  waves  in  the  group  would 
have  a  bed  which  was  over-steep  and  w  hich  would 
efnc;entl\  eject  sand-laden  vortices  high  into  the 
flow  and.  in  the  process,  erode  the  ripple  steepness. 
This  Aortex  generation  from  ripples  which  are  out 
of  equilibrium  (over-steep)  may  allow  a  series  of 
high  waves  to  advect  the  decaying  vortices  high 
into  the. water  column.  Conversely  the  smaller 
waves  would  experience  a  bed  which  was  flatter 
than  at  equilibrium. 


Figures  5b  and  c  (301  3  and  301,11)  show  a 
good  example  of  this  process.  From  ;  =  55  s  to  r  = 
85  s  in  301  11.  virtually  no  resuspension  occurs, 
yet  in  301  3.  where  the  wave-induced  current 
speeds  are  about  the  same,  resuspension  is  occur¬ 
ring  with  everv  wave  and  sand  is  lifted  to  20  cm 
above  the  bed.  We  suggest  that  the  difference  is 
due  to  the  more  energetic  waves  prior  to  /  =  55  s 
in  30!  11  which  have  produced  ripples  of  low 
steepness  (approximately  in  equilibnum  with  the 
higher  waves).  In  the  low  wave  period  which 
follows,  no  vortices  can  be  ejected  because  of  the 
low  steepness  and  the  ripples  then  start  to  build 
in  steepness  (initially  as  rolling  grain  ripples). 
W'hen  the  next  group  of  large  waves  arrives  the 
ripples  are  over-steep  once  again  and  rapid  vortex 
ejection  and  ripple  erosion  occurs. 

8  Conclusions 

Highlighted  by  these  results  is  the  importance 
of  knowing  the  bed  ripple  dimensions  during  sus¬ 
pension  events.  Unfortunately,  photographic  tech¬ 
niques  for  ripple  measurement,  both  video  and 
stills,  are  very  susceptible  to  any  decrease  in  water 
clarity  and  a  tracking  high-frequency  acoustic 
device  is  probably  the  best  method  for  gaining 
suitable  data  during  sand  suspension  events  (Din¬ 
gier  and  Clifton.  1984).  Rapid  profiling  of  the  bed 
in  the  vicinitv  of  the  acoustic  sand  suspension 
measurements  with  a  vertical  resolution  of  a  few 
millimetres  would  enable  many  points  of  specula¬ 
tion  to  be  resolved,  particularly  those  surrounding 
the  variation  in 

The  major  control  on  sand  resuspension  during 
groupy  waves  appears  to  be  that  of  bed  topography 
(ripple  steepness)  being  out  of  equilibrium  with 
the  waves  for  most  of  the  time,  adding  a  further 
complication  to  the  problem  of  modeling  sus¬ 
pended  sand  transport.  This  also  raises  the  ques¬ 
tion  of  the  speed  at  which  ripples  respond  to 
changing  wave  conditions;  the  groupiness  of  sus¬ 
pended  transport  suggests  a  time  scale  of  a  few 
wave  periods  but  if  the  values  truly  show- 
hysteresis  the  time  scale  w-ould  be  very  much 
longer. 

Whereas  due  to  the  variations  in  the  distribu¬ 
tions  of  sand  sizes  with  height  there  must  be 
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considcrabie  uncertainty  in  the  accuracy  of 
translating  the  acoustic  measurements  into  sus¬ 
pended  sand  concentrations,  the  points  presented 
here  should  not  be  substantially  altered  because 
the  assumption  of  a  constant  distribution  should 
leave  the  results  internalK  consistent.  The  measure¬ 
ment  of  the  size  distribution  variations  on  time 
scales  of  an  tndlMduai  wave,  perhaps  through  the 
use  of  rr.ulti-freuuency  acoustics,  and  its  incorpora¬ 
tion  into  the  transport  prohies  must  be  considered 
as  the  next  important  step  in  understanding  and 
predicting  sand  transport  in  the  sea. 
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Abstract 

Field  observations  of  the  suspension  of  sand  near  the  sea  bed  under  shoaled  waves 
in  the  nearshore  region  are  presented  and  discussed  with  reg^d  to  the  influence  of  wave 
groups.  Sand  is  found  to  be  suspended  at  both  the  incident  wave  frequency  and  also  at 
lower  frequencies  corresponding  to  wave  groups.  These  limited  observations  indicate  the 
likelihood  that  wave  groups  result  in  higher  concentrations  of  suspended  sand  than  randomly 
distributed  waves  of  similar  heights. 

Introduction 

The  suspension  of  sand  within  approximately  10  cm  of  the  sea  bed  in  nearshore  regions 
has  been  observed  to  occur  in  an  event-like  manner  over  a  range  of  time  scales  (Bren- 
ninkmeyer,  1976;  Hanes  and  Huntley,  1986;  Sternberg  et  al.,  1989;  Hanes,  1990).  Although 
time  scales  for  near  bed  sand  suspension  comparable  to  the  period  of  the  surface  grav¬ 
ity  waves  were  reported,  the  same  investigators  also  reported  that  suspension  occurred  at 
longer  time  scales  corresponding  to  infragravity  waves  or  wave  groups.  Because  sand  near 
the  sea  bed  in  still  water  would  sink  to  the  bed  in  a  few  seconds,  suspension  at  longer  time 
scales  would  seem  to  require  fluid  motion  at  the  same  time  scales.  In  the  case  of  infra¬ 
gravity  wave  motion  there  can  be  highly  energetic  low  frequency  fluid  motion  which  could 
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cause  sediment  suspension.  Observations  of  this  phenomenon  were  reported  by  Beach  and 
Sternberg  (1987).  In  contrast,  there  is  no  fluid  motion  at  the  wave  group  time  scale  due  to 
the  wave  groups  alone,  so  the  dynamical  links  between  shoaling  wave  groupiness  and  sand 
suspension  are  not  obvious.  One  well  described  and  relevant  hydrodynamic  phenomenon 
related  to  wave  groups  is  the  forcing  of  a  low  frequency  wave  bound  to  the  wave  groups 
(Longuet-Higgins  and  Stewart,  1964).  Shi  and  Larson  (1984)  suggested  that  the  bound 
long  wave  in  combination  with  wave  groups  would  result  in  offshore  sediment  transport. 
This  result  is  based  on  the  reasonable  assumption  that  more  sand  will  be  suspended  during 
intervals  of  larger  waves,  which  will  be  coincident  with  the  offshore  phase  of  the  bound  long 
wave.  However,  no  sediment  transport  models  or  theories  predict  a  group  related  time  scale 
for  suspension  or  an  enhancement  of  suspended  sediment  concentration  during  grouped  se¬ 
quences  of  large  waves,  relative  to  that  which  would  occur  if  the  same  waves  were  randomly 
distributed. 

The  general  physical  linkages  between  the  fluid  forcing  and  sediment  response  are  not 
well  known  under  unsteady  forcing  conditions.  If  suspended  sediment  concentration  were 
instantaneously  proportional  to  the  near-bed  velocity  raised  to  an  even  power,  for  example, 
then  the  low  frequency  envelope  seen  in  the  velocity  record  during  wave  groups  would  be 
reflected  in  the  suspended  sediment  record.  This  partially  explains  the  field  observations 
to  be  presented.  However,  it  will  be  shown  through  spectral  analysis  techniques  that  the 
suspended  sand  concentration  is  actutdly  enhanced  during  wave  groups. 

The  field  measurements  to  be  described  were  obtained  in  the  nearshore  region  at  Cape 
Canaveral,  Florida  with  the  use  of  a  pressure  sensor,  a  2-axis  electromagnetic  fiownneter, 
and  an  acoustic  concentration  meter.  The  flowmeter  w^ls  centered  approximately  15  cm 
above  the  sea  bed.  The  averaging  volume  of  the  flowmeter  was  approximately  1  liter,  and 
the  amplitude  response  W2is  near  unity  for  frequencies  below  1  hz.  Accuracy  in  the  fluid 
speed  measurements  is  better  than  5  cm/s.  The  acoustic  concentration  meter  measured 
the  vertical  profile  of  suspended  sand  concentration.  The  temporal  zuid  spatial  resolution 
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of  the  concentration  measurements  were  0.25  seconds  emd  0.5  cm,  respectively.  Based  upon 
laboratory  calibration,  the  uncertainty  in  the  concentration  measurements  are  believed  to  be 
approximately  10%  for  the  relative  concentration,  and  25%  for  the  absolute  concentration. 
These  uncertainties  do  not  significantly  influence  the  results  to  be  reported  here. 

The  data  were  obtained  over  a  6  hour  period  on  the  evening  of  July  13,  1988,  in  the 
region  just  outside  the  surf  zone  in  water  depth  of  about  2  m.  The  tide  rose  about  25  cm 
during  the  first  3  hours,  and  dropped  about  50  cm  during  the  l2ist  3  hours.  Data  consists 
of  time  series  recorded  at  a  rate  of  approximately  4  hz  for  8.7  minutes  each  hour,  so  the 
duration  of  data  presented  is  approximately  52  minutes.  The  significant  wave  height  at  the 
sensors  varied  about  20%  over  the  duration  of  these  observations.  It  was  approximately 
40  cm  with  a  peak  wave  period  of  8  to  9  sec.  The  meein  cross-shore  current  measured  by  the 
flowmeter  was  7.1  cm/s,  directed  offshore.  The  local  sand  had  a  narrow  size  distribution 
with  mean  sieve  diameter  of  0.18  mm.  The  sea  bed  was  flat,  with  occasioned  shell  fragments 
protruding  up  to  2  or  3  mm  above  the  sand. 

The  data  will  be  described  below  in  both  the  time  and  frequency  domains.  Both  tech¬ 
niques  will  be  used  to  demonstrate  the  correlation  between  near  bed  suspended  sediment 
concentration  and  wave  groups.  The  quantity  is  chosen  to  characterize  the  wave  forcing 
because  of  its  usefulness  in  describing  the  wave  groupiness.  is  also  a  measure  of  the  fluid 
speed,  wave  height,  bottom  shear  stress  magnitude,  energy  density,  etc.  The  dynamiced 
dominance  of  any  one  of  these  quantities  is  not  implied,  and  the  identification  of  specific 
processes  causing  the  correlations  to  be  presented  below  remains  as  future  w’ork. 

A  time  series  of  suspended  sediment  vertically  averaged  over  the  region  0.5  cm  to  4.5 
cm  above  the  sea  bed  (NBC)  is  shown  in  Figure  1  along  with  U^,  where  U  is  the  cross¬ 
shore  component  of  the  fluid  velocity.  Figure  1,  which  shows  the  first  of  the  six  data  runs, 
illustrates  the  previously  mentioned  observation  that  suspension  occurs  at  a  time  scale  near 
that  of  the  incident  waves,  and  also  at  longer  time  scales.  The  longer  time  scales  are  further 
illustrated  in  Figure  2,  which  shows  lowpass  (30  sec)  filtered  time  series  of  both  and 
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NBC  for  the  entire  six  records,  plotted  in  sequence  with  a  sUding  ordinate.  The  visual 
correlation  at  long  time  scales  is  striking.  Correlations  to  be  presented  below  are  calculated 
as  follows.  NBC  eind  17*  for  each  data  run  are  lowpass  Sltered  and  cross  correlated.  The 
maximum  correlations  are  then  averaged.  The  results  are  presented  in  tabular  form  below: 


Lowpass  Filter 

Correlation 

95%  Confidence  Limits 

None 

.30 

.29  to  .31 

30  second 

.71 

.59  to  .79 

100  second 

.67 

.41  to  .83 

The  confidence  limits  are  estimated  zissuming  normal  distributions  with  12228,  102,  and 
31  degrees  of  freedom.  These  correlations  are  remarkably  high;  they  strongly  support  the 
hypothesis  that  wave  groupiness  is  causally  linked  to  near  bed  seuid  suspension. 

Figure  3  presents  the  power  spectrum  of  NBC  and  17*  along  with  the  95%  confidence 
intervals.  Individual  spectra  were  computed  using  1024  point  fft’s  with  no  overlap,  and 
then  averaged  together.  In  contrast  to  the  spectrum  of  [/  (not  shown)  which  has  one  to 
two  orders  of  magnitude  less  variance  at  the  group  periods  compared  to  the  incident  band, 
the  spectra  of  17*  and  NBC  show  considerable  variance  at  low  frequencies.  As  alluded  to 
previously,  one  could  argue  that  the  suspended  sediment  concentration  is  instantzmeously 
proportional  to  17*,  and  that  suspension  at  group  time  scales  simply  reflects  17*  at  group 
time  scales.  If  the  NBC  were  perfectly  linearly  related  to  17*,  then  the  spectral  ratio  between 
NBC  and  17*  would  be  constant.  Rather,  Figure  4  shows  that  this  ratio  increases  as  much 
as  a  factor  of  2  to  3  in  the  50  second  to  100  second  range  of  periods. 

Figure  5  shows  the  spectral  coherence  (computed  as  before)  between  17*  and  NBC,  where 
[//is  the  high  pass  filtered  U  with  a  cutoff  of  20  seconds.  The  purpose  of  examining  17*  is  to 
remove  any  low  frequency  fluid  motion  which  may  be  unrelated  to  wave  groups.  The  high 
coherence  between  17*  and  NBC  over  a  broad  range  of  low  frequencies  clearly  demonstrates 
the  effect  of  wave  groupiness,  isolated  from  direct  (non-group  induced)  fluid  forcing  at  low 
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frequencies.  Coherences  above  the  dashed  line  indicate  with  95%  confidence  that  NBC  and 
Uj  are  not  linearly  independent. 

In  conclusion  we  have  identified  two  features  of  a  set  of  field  observations  of  near  bed  sand 
suspension  under  shoaled  wave,  flat  bed  conditions.  Firstly,  there  is  a  significant  correlation 
between  near  bed  concentration  and  cross-shore  fluid  speed  squared  over  a  range  of  time 
scales,  and  secondly,  the  concentration  is  enhanced  at  time  scales  corresponding  to  groups 
of  incident  waves. 
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Figure  3:  Spectra  of  NBC  and  U^. 


BC  to  t/2. 


